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I ^ 1 1 Abstract 

I [ This paper is a continuation of [Brlvrj and |Ivrll llvrSl IIvr6l lIvrTl llvrSj . I consider 

^ ' two-dimensional Schrodinger operator with degenerating magnetic field and in the 

■ generic situation I derive spectral asymptotics as /i — > +0 and fi — > +oo where h 
, and /x are Planck and coupling parameters respectively. The remainder estimate 
I is 0{fi~^ h~^) which is between 0{fi~^h~^) valid as magnetic field non-degenerates 

■ and 0{h~^) valid as magnetic field is identically 0. As fx is close to its maximal 
, reasonable value 0{h~^) the principal part contains correction terms associated with 
I short periodic trajectories of the corresponding classical dynamics. 

■ 

^. ■ Introduction 

>< ■ 0.0 Preface 

H : 

. - - . We consider spectral asymptotics of 

(0.1) ^ = P,9''ix)P, - V) , P, = D, - ^V, 

where g^'', Vj, V are smooth real- valued functions of x G and (g^^) is positive-definite 
matrix, < /i ^ 1 is a Planck parameter and /i ^ 1 is a coupling parameter. We assume 
that y4 is a self-adjoint operator. 

*Work was partially supported by NSERC grant OGP0138277. 
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In contrast to my recent papers |Ivr6[ IIvr7| llvrSj I assume that all the coefficients are 
very smooth and in contrast to |Ivr7j I consider only two-dimensional case here. However 
degeneration of the magnetic field makes result much more interesting and difficult than I 
expected when I started this work. 

I hope to investigate similar even-dimensional case in the future; odd-dimensional results 
are much easier and of no special interest because only in the even- dimensional full-rank 
case magnetic field can improve the remainder estimate. 



0.1 Assumptions and notations 

Let g = det{g^'')-\ F12 = d^,V2-d^^Vi and / = Fug'^V-K Note that both F^g'^ and / 
are coordinate independent. In |lvr6j and in earlier papers I assumed that F disjoint from 
but I will not assume this anymore; however I assume that degenerations are generic, 
namely 

(0.2) V > eo, 

(0.3) IF12I + IVF12I > eo, 

(condition ()0.2j) will be dropped in the very end by rescaling arguments). Then S = {x : 
-^"12 (a;) = 0} is a smooth curve. One can introduce a coordinate xi as a Riemannian distance 
from X to S; then g^^ = 1. One can always change X2 to keep g^^ = and then 

(0.4) A = ^ + p^g~^p^ _ , P, = D, - /iV,-. 

Furthermore, without any loss of the generality one can assume locally that 
(0.5) 1^1 = 0, \/2x|xir 

with z/ = 2; one can achieve this by the gauge transformation. Actually instead I consider 
operator ()0.4j) with positive integer u > 2 and in the classical dynamics I even take any 
real u > 2. 

What I am interested is an asymptotics as h ^ +0, fi +00 of 
(0.6) /e(..x.O).)(x)dx 

where e{x, y, r) is the Schwartz kernel of the spectral projector E{t) of operator A. I will 
prove later that ()0.6|) will be 0{fi~^) as /i > Coh'" and therefore I will assume that 

(0.7) /i < CqH-". 
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The natural answer coming from the non-degenerate case is J S^^{x, 0)ip{x) dx where 
(0.8) = ^^9{2t + V - {2n + l)fxhF)fxh-'Fy^ 

n>0 

(with 6{t) = 0,lasr<0,r>0 respectively) which implies that there are actually two 
different cases 

(0.9) 1 < < h-\ 

and 

(0.10) h^^ <fi< Coh^" 

when the answer will be of magnitude and (nh)^^^^'^~^^h~'^ respectively; in the latter 
case the contribution to the answer will be given by the strip 

(0.11) Z = {x: <7i=Co(/i/i)-^}. 

The standard rescaling procedure applied to results of |Ivr6j implies 

(0.12) n=\ J(^e{x,x,0)-S^'^{x,0)^tfj{x)dx\<Ch-^ 

but an aim of this paper is much better estimate, up to 0{n~^h~^). 
I am going to improve this remainder estimate. 

0.2 Results 

Theorem 0.1. Let conditions (j^.jjj) . i\U.!^ and ()^.5|) be fulfilled and let ip he supported in 
the small enough neighborhood of {xi = 0}. Then 

(i) As 1 < fi < = Ch^"^^ estimate TZ < Cfi^^^^h^^ holds; under some nondegeneracy 
assumptions one can push ^\ up; 

(a) As h~^/'^^ < fJ'2 = Ch"^ estimate 

(0.13) 7^* = I j (e(x,x,0)-£:^W(x,0))^(x)dx-y" (^2, 0)^^(0, X2, 0) rfxal < 
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holds with arbitrarily small exponent S > 0. Here and below correction term is defined 
by l\3.4^ in the terms of the eigenvalue counting function no of the axillary 1-dimensional 
Schrodinger operator with semiclassical parameter h = jjL^^'^h. 

(Hi) Under non- degeneracy condition |(2.105)^ for W{x2) = 1^(0, X2) (see below) as ^ < 
eh~^ with small enough constant e = e(eo,rra) (where eo is a constant in |(2.105)^[ ) and 
under non- degeneracy condition ( 4-151 as eh~^ < < Ch~^ one can skip the last term in 



the right-hand expression of (\U.1'J^ : 
(iv) As fi > Ch~^ |e(x, X, 0)1 < C"/i~* with arbitrarily large exponent s. 

Remark 0.2. With some error making estimate sometimes not that sharp one can replace 
no by Bohr- Sommerf eld approximation for it, providing formulae ()3.52|1 . |(3.52)*[ |(3.52)** 

for ^r"' 



•'corr 



0.3 Plan 

I start from Classical Dynamics which is useful both for understanding and proofs of our 
main results. In Section 2 I consider corresponding quantum (semiclassical) dynamics 
and the remainder estimates. It appears that there will be forbidden zone > 71 = 

C{iihY/^''~^^} outer zone Z^nt = {70 = Cfj.~^^'^ < \xi\ < 71} and inner zone Zi^n = 

{l^^il < 7o} divided in non-periodic zone and periodic zone with the analysis in the latter 
one the most difficult and interesting. The width (in the phase space) of the periodic zone 
is defined by the uncertainty principle and it increases as fi grows. 

In Section 3 I instead of implicit main part of asymptotics derived in Section 2 provide 
much more explicit answer S^"^ + . 

All the analysis in sections 2,3 is under assumption fi < h~^^^ while section 4 is devoted 
to the similar analysis in the case of superstrong magnetic field h~^^^ < /i < ^h~^ which is in 
some sense easier and in some sense more difficult than the previous one because h is not very 
small anymore. Also in Section 4 I consider ultra-strong magnetic field eh~^ < /i < Ch~^ 
when h is not even small anymore; as I mentioned case fi > Ch~" is trivial. 

Finally, Section 5 is devoted to generalization for vanishing V and discussion of other 
generalizations. 

Only as Ch^^ < < Ch^"; otherwise 71 is artificially set to e. 
Here and below coordinate xi is defined so that E = {xi = 0}. 
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1 Classical dynamics 

In this section I consider classical dynamics on energy level described by Hamiltonian 
(1.1) « = 9"'ix)PjPk - V) , p, = - f^V, 

where all conditions of subsection 0.1 are assumed to be fulfilled. 

1.1 Pilot-model 

Let us consider first a model 

(1-2) « = ^(e? + (6-^kir)^-i) 

with V >2, corresponding to yU = 1, V = 1. 

Proposition 1.1. For symbol () j.jjj) 

(i) Along trajectories ^2 = const, and 

(a) if ^2 > 1 on energy levels] these trajectories oscillate either between Xi = hi andxi = 62 
or between xi = —hi and x\ = —62 with hi = ((^2 — 1)^)^^''? ^2 = ((G + 1)^)^^'^/ 

r^y' ^/^2 = 1 these trajectories vary between hi = and 62 = 1 or between —1 and 0; there 
is also exceptional trajectory xi = = 0, X2 = t + const; 

(c) if ~1 <^2 < 1 these trajectories oscillate between Xi = hi = —62 (in'd Xi = 62/ 

(d) if C,2 = —1 energy level degenerates into {xi = ^1 = 0} and this trajectory is 
Xi = = 0, X2 = i2t + const; 

(e) energy level is empty as ^2 < —1/ 

(ii) Along each of these trajectories (xi,.^i) are periodic with period 

(1-3) T{i2) = 2 r —==L=== dy 



lb, Vl - (A; - lylVz/) 
and X2{t) = v{^2)t + 5:2 (^) where v = T^^I , 



b2 



k - |y|7z/ 



Jbi ^Jl-[k-\y\''|vY 
and X2{t) is T{^2)-P^™dic function. 
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Proof. Proof is obvious consequence of the fact that evolution in is described by 1- 

dimensional Hamiltonian p.2j) while ^2 = k = const. Potential W = W{k) = (/c — — 1 
has two wells as /c > and one well as /c < 0; as A; = it has a fiat bottom. Note that the 
bottoms of the potential are in ±6(/c), b{k) = [kuY^'^ . 




(a) k = 1.5, two wells (b) fc = 1, two touching wells (c) k ~ 0.9, one well 

Figure 1: Graphs of W{k), 1-st pilot-model 

Further, as /c > 1 we have W{Q) > and thus the particle lives in one well oscillating 
either between bi(k) = ((k — l)vY^^ and h2{k) = {{k + l)iyy^'^ or between —h2{k) and 
—hi{k). As < 1 we have either one-well potential or W{Q) < so particle oscillates 
between hi{k) = —h2{k) and 62 (fc)- As = 1 particle moves between hi{k) = (as t — * ±00) 
and ±62 (^) (as long as z/ > 2). □ 

Proposition 1.2. (i) v{k) > and both T{k) and v{k) are decreasing functions of k > 1. 
(11) The width of the wells is ~ 26(1"'^)/'^ as k ^ +00; T{k) ~ 2Txh^^'' , v{k) ~ \{u - l)h''' 
ask^ +00, h{k) = ikpf/''; 

(Hi) T ~ 2 1 log l/c — 1| I (for v = 2), T const|/c — (for v > 2) and ~ 1 as /c— ^1+0. 
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Proof. Changing z = y'^/ u — k we arrive to 

(1.5) T =2 ((fc + z)u)T^~\l - dz, 



(1.6) ^"~^y ^(('^ + ^)^)''"'(l-^')~^^'^ = 

' z(^{{k - z)u)^-^ - {{k + z)u)^-'^ (1 - z^)-! dz, 

which imphes that I > as k > 1. Asymptotics T ~ 2-K{kv)^^^ , I ~ 7r(z/ — l)(A;z/)^~^ as 
— > +00 follow from ()1.5j) . ()l.fi|l as well. 
Also ^^,^Ms imply that 

(1.7) ^ = -2{v-l) j'^ {{k + z)v) '-'(1 - ^2)-^ dz, 

(1.8) |^ = - 2(1.-1) j\{{k + z)uy-\l - z^)--^ dz, 

(1.9) 1^ = - 1.(1. - l)T-2 /■\^_^^)2((^ + ^)^)^-2^(fc + ^^)j,)^-2^ 



1 >/-i 

2 



{l-z^)-^dz, 

implying monotonicity properties. 

Asymptotics as — ^ 1 + follow directly from ()1.3p . (11.41) . □ 

Remark 1.3. In the strong magnetic field the drift is described by \Cf^^ = {dx2f, —dxif) = 
(0, {i' — l)xi'^) (where | comes from this factor in the front of Hamiltonian) and this matches 
to asymptotics of f as ^2 — ^ +00. 

As — 1 < A; < 1 one can rewrite 
(1.10) T = Aj'j{l-{k- y^lvf) dy = 

A{k + (^2-{k + 1)(1 - y'')Y" dy = 

4{k + l)^-'u^ [ {l-y''y'^^{u,y'')dy 
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and 



(1.11) I = Aj^ ^k-y''/v){l-{k-y''/ufY' dy 



A{{k + l)u)"- [ {l-y-r-^^^{uy)dy 
Jo 

with u = k/{k + 1) e {—oo, |), z) = {1 — 2u + z)~^^ V^i(m, z) = {u — z){l — 2u + z)~^. 

Let Ti{u), Ii{u) denote integrals in the right-hand expressions of (jl.lUj) . (Il.llj) . Then 
g > 0, = (1 - > because |^ = (1 - 2n + ;.)-§, ^ = (1 - n)(l - 2n + z)-l. 

Also 

due to Cauchy inequality and therefore 

I = - + > + - '''^') = ^ i"'^ - - + > 0. 

Proposition 1.4. ("zj T/iere exists k* = k^ e (0, 1) snc/i that v{C,2){C,2 — k*)~^ > and it is 
both bounded and disjoint from as k E { — 1 + e,l) ; 

(a) v{k) is monotone increasing at (—1, 1); 

(Hi) T{k) ~ consi{k + 1)^ and v{k) ~ — 1 as k ^ —\ + 0; 

(iv)T ^ 4|log|fc-l|| (foru = 2), T ~ 2const|A;-l|^-^ (foru > 2) andv ^ 1 as k ^ 1 + 0. 

Proof, (ii) is already proven; (iii),(iv) are easy and (i) follows from (ii)-(iii) and from f (0) < 
0. □ 

Remark 1.5. (i) T{k) is monotone decreasing function of k for z/ = 2; not sure if this is 
the case for larger u; 

(ii) ^ > iff |xi| < {kv)^] in particular ^ < on the whole trajectory iff — 1 < < 0; 

(iii) One can prove easily that || ~ =pconst|A; — 1|^^~^| log \k — 1||~2<5^2 as A; ^ 1 ± 0; 

(iv) In particular, X2{t) is periodic iff C,2 = k*; 

(v) Since 

^ {u-l)'{l-yTkl-2u + yTh%-logy)dy>0 



dv 

we conclude that /c* monotonically decreases; one can see easily that A;* — > +0 as z/ — > +oo; 
(vi) Maple experiments show that k2 ~ 0.65 (but not 2/3). 
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The following series of figures El show (a;i, X2)-trajectories at level for different values 
of .^2 = ^ as z/ = 2. In particular figure p(g)1 shows periodic curve. 

Proposition 1.6. There exist functions Z = Z{xi,C,i), a = a{xi,C,) and [3 = P{xi,^2), 
such that 

(1.12) {a, X2-Z} = ^2- {a, Z} = ^(6 - k*) + (3a{xu 0- 
Further, Z is odd and a, (3 are even with respect to each of Xi, ^i. 

Proof. Consider now periodic curve as .^2 = k* with the "center" at 0. It is described by 
equations in (x, .^)-phase space 

(1.13) G = k\ a(xi, = 0, X2 = Z(xi, ei) 

with smooth odd with respect to each xi, ^1 function Z. We cannot express it via Xi or ^1 
alone. Then 

{a, X2 - Z} = ^2 - {a, Z} = as ^2 = k* , a{xi, ^) = 0. 

Note that da and d^2 are linearly independent unless .^1 = and (.^2 ~ \xi\'^/i')xi = which 
for a = means that ^2 = ±1- Therefore we arrive to ()1.12|1 . □ 

Remark 1.7. For even integer u all these functions are analytic. 

1.2 Another pilot-model 

Consider now another model 

(1.14) a = 1(^1^ + (6 - i|xir sign(x))' " l) 

with u > 2, corresponding to fi = 1, V = 1. Again we are looking at energy level 0. 
Consider again the graph of potential. 

In the similar way movement in (xi,^i) is described by Hamiltonian with one-well 
potential and if/c>lor/i;<— 1 this particle moves only in {xi > 0} or {xi < 0} 
respectively and in the former case the evolution coincides with given in the previous 
subsection while in the latter it one needs to replace k by —k, take a mirror-symmetric 
picture but now the directions "up" and "down" are switched while direction of rotation 
does not change. So outer zone is described exactly as we did (see figures ^-d) but in the 
inner zone for — 1 < A; < 1 we get a kind of combined picture (see figure E}. 
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(c) k ~ 1.1, w w 1.13 (d) k ~ 1, w 1.13, first critical case 

Figure 2: Drift up; and rotation clockwise, each figure has its mirror with respect to a;2-axis 
with drift up and rotation counter-clockwise; w = \/k + 1 — \/k — 1 is the width; 
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Figure 2: (continued) As k decreases drift up changes to drift down; on |2 (g)| movement is 
periodic; w = 2\/k + 1 is the width; 
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(k) k = -0.5, w w 1.41 (1) k = -0.9, w w 0.63 

Figure 2: (continued) As k decays movement is straighten; drift is down. 
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(a) k = 1.5, well in xi > (b) fc = 1, well touches xi = (c) fc = 0.9, well contains 

XI = 

Figure 3: Graphs of W{k), 2-nd pilot-model 

Then 

T{k) = 1(T+(A:) +T+(-A:)), I{k) = i(/+(A:) - 

where Tj^{k), I+{k) are exactly what we introduced for Hamiltonian ()1.2|) . However due 
to monotonicity of ||/|'' sign(?/) we can apply idea of proposition 11.21 introducing z = 
|y|'^sign(?/) — k we arrive to 

1 
1 



;i.i5) T{k)= 2j \{k + z)i^\"-~\i-z^yUz, 

[1.16) I{k) = -2 j' z\{k + z)u\^-\l - z^)-^ dz 



-1 

2 



z[\{k - z)v\^-^ - \[k + z)i^|^~')(l - z')-^ dz 



and obviously eQl{k)/k < cq as 7^ 0. 
Thus we arrive to 



Proposition 1.8. Consider on level trajectories of Hamiltonian [ 1.14 ) with $,2 = k & 
(—1,1). Along them 

(i) Variables (a;i,^i) are T(k)-periodic and Xi oscillates between hi = —((1 — k)uY/'' and 

62 = ((1 + 
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(c) k = 0.1 (d) = 

Figure 4: Drift is up and movement clockwise. As k decays from 1 to drift up slows 
down; |4(d)| corresponds to periodic case; we skip k > 1 which repeat those at figure Efe^-d, 
A; < — 1 which will be mirror-symmetric to them with movement down and also clockwise 
and — 1 < A; < as symmetric to given there. 
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(a) X2{t) = X2{t) + v{k)t with T-periodic X2{t) and v{k) = I{k)/T{k) ; 
(in) eoI{k)/k < cq as k ^ 0; 

(iv) Therefore the only periodic trajectory is with k = 0; 

(v) Proposition ] l.b] holds with k* = but Z is even and a,/3 are odd with respect to ^i; 

(vi) For odd integer u > 3, Z , a, (3 are analytic. 

1.3 Pilot-models perturbed 

The periodic trajectories of our pilot-models are very fragile and one can destroy them 
easily. Consider heuristically the same Hamiltonians with V = ax\ with a > instead of 
V = \. Then ^ = \dx.y = \oi and thus ^2 = cq + \at. Then the averaged movement 
along 1/ = 0:2 is described by 

'Ly = V^v{^2{t)V-^) ~ {ay)h!{^at{ay)-^) j^z ~ v{tz-^) 

with z = l^yjat. This "equation" has a solution z = (3^^t as ±t > as /3±v{f3±) = 1 and 
±P± > 0. 

In the second pilot-model this equation has solutions with — = > 1 and therefore 
on them "in average" Xi ~ {u^2V~^)~y^ ~ py~ ■ In the first pilot-model /?_ = —1. These 
observations explain the evolutions in vertical direction on figures Eli^-b. 

1.4 General case. Remarks 

We are going to consider 

(1.17) a = ^ (e? + (j{xY{^2 - /i0(x)^?,(xi)/z/)' - V(x)) 
with coefficients 

(1.18) a>eo, 0>eo, V > eo 

and all of them are smooth and with = \xi\'^ or g^, = signxi, z/ > 2. 
Let 

(1.19) Pi = 6, P2 = ^2- fJ'(i){x)g^{xi)/iy; 
then 

(1.20) {pi,P2} = -l^d^i{<P{x)Qu{xi)/i^) = -iJ.(p{x)g,,^2{xi)xi{l + 0{xi)). 
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(a) First pilot-model perturbed (b) Second pilot-model perturbed 

Figure 5: Both pilot-models are perturbed by a weak electric field directed up (V^ = ax2) 



We want to consider trajectories on level {a = 0} (and thus on any level {a = r}) and 
thus we can multiply a by any smooth factor disjoint from 0; then we can assume that 

(1.21) a\^ = ^\^ = l. 

Really, starting from (jl.lTj) we can multiply a by smooth function ci;^ and redefine Xi 
according to the new metrics; then we redefine X2 to get rid of the mixed term. Then 
fll.l7|) will still hold but with V i-^ u^V, a\-s i-^ (crcij^)!^, 0|s ^ {(j)uj'-')\^ and picking up 

u = (cr0)~^/*^'^+^) we arrive to the case (cr</))|s = 1 with 

(1.22) V^V* = (cr0)~2/('^+^V. 

Finally, replacing X2 by p{x2) we can achieve a|s = 1 and thus 1)1.211) . 



1.5 General case. Outer zone 

Proposition 1.9. For general Hamiltonian ( |i.i7| ) in B{0, 1) satisfying 1^1.18^ 
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(i) There exists b = b{x,pi,p2), such that as \a\ < and |/| > 7 = C/i ^^'^ 

(1.23) Co' <J<Co, 

(1.24) \{a,b}\ < Cofi^'Ur'; 

(a) Trajectories starting at t = from {x,C,) G 5(0, |) fl {x,±f{x) = 7 > 7} for time 
\t\ <Ti = eo/i/(x)^ remain in B{0, 1) and along them e < ±f{x{t))/'-f < C; 

(Hi) Trajectories starting at t = from {x,C,) G i?(0, ^) fl {x, |/(a;)| < 7 = 7} for time 
\t\ < Ti = e remain in -8(0, 1) and along them \ f{x{t))\ < C17. 

Proof. Let us consider trajectory starting at (x,^) with x G -8(0, |), a < c, 7 = |a;i| > 7 = 
Co/i^^/*^. Consider this trajectory as long as it is contained in -8(0, 1) but not further than 

Along this trajectory 

(1.25) ^ = {a, 6} = a\x)p2(l).,ix) x fig,ixi)/u + 0(1). 



Then for X = ^ - (]{x)p 



dX 1 

(1.26) — = {a,X} = a^P2^Qu{x,){-4>.,{x) - <P{x)(3x^^) + 0(1) 

(as a = 0) and therefore = 0(1) for 

(1.27) X = ^2-P{x)pi, /?= = ixi9,,log|0(x)|. 

Therefore |X(a;) — X[x)\ < CT < Cenj'^, which together with inequality \p2\ < c imply 
that |.^2 — ^2! < Cefi'-y"] for small enough constant e > we also get eifi'j'^ < \C,2\ < cifi^y" 
(I remind that fj.'j'^ > C) and thus e^'^'j < |a;i| < C17. 

So, trajectory remains in the strip described above and since the drift speed according 
to |Ivr6j does not exceed Co{fi'y'^)~' there we are insured that trajectory remains in -8(0, 1) 
for time T. This proves (ii). 

On the other hand, if trajectory starts in {x,^) with < 7, it remains in -8(0, 1) for 
time T = e since the speed does not exceed Cq and due to above arguments it cannot get 
into zone {xi > O7}. This proves (iii). 

To prove (i) consider b = (/i^^X)^/*^; then as \xi\ > 7 obviously 6 x |xi| and furthermore 
|{a, b}\ < 0(/i-^|6|)^/^|^2r^ < cjj-^lxil^-" which is exactly (fT^ . □ 

3)Then Ti x 1 as \xi\ x 7 and Ti x ^ as \xi\ x 1. 
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Obviously = <j'^P2 and for F = X2 — /i ^ ()^i,{xi)xil3'pi we have 
rIY 

So, let /5' = then 

^ = - l)/3V-if,_,(xi)(a - + ^(z/ - 1)(3'^^~'q..{xM - a'pl) + 0{^^^'\x,\^-n 
and redefining 

(1.28) F = X2 - ^^~^Q~v{xl)xlj3'pl - fi'"^ g^2u{xi)xiPiP2 
we arrive to 

HY 1 

(1.29) — = --{u - l)(3'fi~'g.,{x,) + 0{fi~'\x,r''' + ^I~'\x,\'"'') , 
and we arrive to 

Proposition 1.10. Along trajectories in outer zone {e > > 7} 

(1.30) eo{firr'\t\ < ±goix){Y{t)-Y{0)) < C{fiYr'\t\ 
asO<±t< efiY- 

1.6 General case. Inner zone 

Now we need to consider trajectories lying in inner zone Z^^^ = {x : \xi\ < C7}. Here 
function X is still defined but Y is not and according to subsection 1.1 the important role 
is played by X{x,p)V{0,X2)~^. 

Proposition 1.11. Let conditions i\l.lS^ . |^2(0)| < C be fulfilled and let condition 

(1.31) '^(6-^*1^1^) > ei, '^ = ±1 
be satisfied as t = 0. Then 

(i) For \t\<T = e this condition is satisfied (with ei replaced by €i/2); 
(a) Moreover 

(1.32) e2\t\ < ±<;{x2{t) - X2i0)) < C2\t\ as c^<±t<T 
and 

(1.33) e2\t\ < \x{t) - x{0)\ + \p{t) - p{0)\ < C2\t\ as |t| < 07. 
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Proof. Proof of (i) is obvious; ()1.H2|) and p.HHj) then follow from comparison of the general 
and pilot-model system. □ 

So there are no T-periodic trajectories with < T < e unless 

(1.34) P=\^2-k*V^\<ei. 
Note that now 

(1.35) f-^'^^ + OW 
which immediately yields 

Proposition 1.12. Let conditions 1^1.1^ . 1^1.21^ . I\l.,y4l be fulfilled and 

(1.36) <;d^^V > eo with ^ = ±1 
Then trajectories on level satisfy 

(1.37) ±^(6W - 6(0)) > e\t\ as < ±T < e. 

In what follows we will consider more precisely trajectories satisfying ()1.34|) and also 

We still want to consider such trajectories as long as they are in a certain vicinity of the 
original point x(0). For the pilot-models this vicinity is B{0, 1) and T = ep~^ as 1^ — = p 
but for other Hamiltonians both vicinity and time are smaller. 

Proposition 1.13. For any of two pilot-models with coupling constant p> 1 and potential 
y = 1 ± a, < C < e 

(i) Trajectory starting at (x,^) satisfying J^ ) remains in -8(0, 1) and satisfies ( 1.34 ^ with 
2ei as long |t| < T = emin((^~2, p~^^; 
(a) Along it 

(1.38) \x2{t) - X2(0)| < + P\t\) + C7, Ut) - 6(0) = ±^Ct 
and in at least one time direction 

(1.39) \x2{t) - X2(0)| > e(C|tp + p\t\) - Cr, 
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Proposition 1.14. For general Hamiltonian j?D , satisfying ()i.i<^ . (ji.j^ij) anc? 
(1-40) \d.,V\j:\<C 

(i) Statement (i) of proposition f7T73l holds as 

(1.41) p > 7, C > 7 

anc? — k*V^ \ < p; 

(ii) Assuming ^l.^l ), if 



|e2(t)- 6(0)1 >eC|t|; 



(1.42) 

then 
(1.43) 

(Hi) Assuming ), 

(1.44) |6(0)-FV(a;(0))^| >ei/7 
i/ien m at least one time direction 

(1.45) \x2{t)-X2m>ep\t\. 

Proof. Proofs are obvious. □ 
1.7 Hamiltonian maps 

Now we need to analyze more precisely what happens in the periodic zone Z^^v = {x-, \xi\ < 
cii — k*V'i\ < p} where p will be specified in the next section. Let us consider pilot- 
models first. 

If on each energy level r, |r| < e all the trajectory were periodic with period T = T{t) 
then replacing Hamiltonian a by g[a) with gij) = J^T{T)dT we would get 1-periodic 
Hamiltonian fiow. 

However it is not the case: on each energy level r the only periodic trajectory is one 
with 6 = 6(r) = k*{l + 2ar)i and period is T(r) = T*(l + 2r)(''-i)/2'^. Still defining ^(r) 
accordingly let us consider the Hamiltonian g{a). 
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Proposition 1.15. For a pilot-model symbol a given by (I i.jjjl or (\1.14l 

(1.46) e^*-^f('') = 

(1.47) ^(-)-^(l + 2r)^-^' 

(1.48) 6 = 6(6,^) = (1 + 2a)" f ^^-k*Yuj( 



2a)5 / V(l + 2a)^ 
where T* is an elementary period and 



:i.49) o;(r) = i9,/(A:)|,, 



k=k* 



Proof. As a is a pilot-model Hamiltonian then taking C,2 = k*, a = we would conclude 
that the Hamiltonian trajectories are periodic with period T*. As a = t disjoint from — | 
we have both period T(r) and K*{t) defined by 

(1.50) T(r) =T*- (l + 2r)(^"^)/(2^), i^*(r) = F ■ (1 + 2r)5. 

To get period which is independent on energy level one must replace a by g{a) with g{r) = 
T{t) dr: then e^*^s{a) = / as ^2 = K*{a)] period is T* rather than 1 since we define g{a) 
without factor T*. 

As ^2 7^ K*{a) periodicity is broken. Since Ha and H^^ near {a = 0} are linearly 
independent we conclude that = J + 2(^2 — K* {o,))uJi{xi, ^i, ^2) ■ However since it 

is a symplectic map it must be of the form e^'' with 

6= (e2-i^'*(a))2^2(xi,ei,6)+A(xi,ei). 

We know that as .^2 = K*{a) we have e^'' = e^^ = I and since A is "small" we conclude 
that A = const (and thus we can take it 0). 

Since e^'' commutes with we conclude that u = 00(^2,0,) and due to homogeneity 
properties UJ2 = Cij(^2(l + 2a)~2^(l + 2a)'^~^. Then X2-shift as a = and ^2 ~ k* is 
2uj{k*){^2 - k*) +0{{^2 - k*f) which implies flT3n|l . □ 

If we consider pilot-model but with /i 7^ 1 one needs to replace ()1.46|) by 

(1.51) e^^*^(") = e^"\ 

Consider now trajectories residing in 5(0, C07) and corresponding to these Hamiltonians 
but with potential V = 1 + C,X2- One can see easily that ()1.51|) will remain true with 
b = b — ^T*W + b', W = V{0, X2) and 6'satisfying ()1.53|) below. Then for the general 
Hamiltonians one can prove 
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Proposition 1.16. For general symbol j?D with a, (f) satisfying !^1.21^ and with V , V{0) = 
1, satisfying in Cq'J -vicinity of 

(1.52) \d,,V\ < C 

with C > C7, equality () j.5ip holds with symbol b = b — ^T*W{x2) + b' where b = b{C,2, cl) is 
defined by ( |i.^<§D , and 

(1.53) \d:b'\ < C(C(C + P) + 7)f l«l < 1- 



2 Quantum dynamics 

2.1 Preliminary notes. Forbidden zone 

Starting from this section 1/ > 2 is an integer and magnetic potential is (0, 0(x)f)y(xi)/z/) 
where Qy{xi) = x\. Thus we have different pilot-models for even and odd v. 

In this subsection we consider /i > and forbidden zone is defined by p.ljl . Also we 
consider ,^2-localization which is useful when /x is close to h~'^ . 

Proposition 2.1. Let fi > and 

(2.1) l^hf-V>e 
in B{x, 7) with 

(2.2) 7 > max(/;,, 7) , 7 = C/i~~. 

(2.3) |e(x,x,0)| < C7"^(/i~^/i7~'^"^)' «n 5(x,7). 
Proof is just by rescaling. In particular, as 

(2.4) h-^ <fi< Coh-" 

the total contribution to J e{x, x, 0)^/'(x) dx of the forbidden zone defined by ()2.H) does not 
exceed C'ji^{fih'^y. Thus we need to consider only allowed zone jjihf — V < e and there 

(2.5) 7<7i = Co(/i/^)"^. 

This result is not very useful as /i > h~'^~^^. However using logarithmic uncertainty 
principles and related arguments one can prove easily 
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Proposition 2.2. Let conditions . be fulfilled and 

(2.6) 7 > max(C/i| log/i|,7). 
wzi/i C = Cs- Then 

(2.7) |i^t^h-irXT(t)V^M| < CTh' Vr < e VT > log/i| 
and therefore e{x,x,0) < Ch^ in B{x,'~f). 

In particular, ()2.6p is fulfilled automatically as > 71 provided 

(2.8) /i < e(/i|log/i|)-^ 



while the super-strong magnetic field case when fj2.8|) is violated requires some special 
consideration based on the analysis operators with operator-valued symbols. We need 
axillary 

Proposition 2.3. Let A = A„(2) be eigenvalues of one- dimensional Schrddinger operator 

(2.9) Si{z)'' = Dl + {z-fig,{xi)u/uy 

in L^{M.); here we assume only that 2 < z/ G M and either g,^{xi) = \xi\'^ or g^i^Xi) = 
■ sign(xi). Then 

(i) \n{z) > 0; 

(a) Assume that either z >0 or = sign(xi). Then for n < e|2|^ 

(2.10) Xn-lzl'^n, 

d 

(2.11) TT-^n ■ sign(2:) x \z\~~n; 

oz 

(Hi) Assume that z < and g^, = \xi\'^ . Then 

(2.12) Xn-z'^+n^. 



Proof. Proof of (i) is obvious; (ii),(iii) easily follow from the semiclassical character of the 
spectrum as |2;| + n ^ 1, n ^ 1 respectively. Easy details are left to the reader. □ 

Furthermore, using operators with operator-valued symbols one can prove easily 
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Proposition 2.4. Let fi > Ch"". Then 

(2.13) \e{x,x,T)\ <Cfi~' for r < e(/i/i^)^. 

So, in what follows can assume that fi < Ch~^ . In what follows u = u{x,y,t) is the 
propagator, i.e. the Schwartz kernel of operator e*^ 

Proposition 2.5. Lei /i < Ch-" and let ip E C~(M), ip{t) = 1 ast>2, ip{t) = as t < 1. 
Then 

(i) The following inequality holds 

(2.14) \Ft^h-^,XT{t)^{±e{nh'')^^hD^^)%l)u\<C{T+l)h' Vr < e VT > log /i| 
and therefore 

(2.15) |(/?(±e(/i/i'')^/iL)2)e(x,y,0)| <Ch'] 

here and below x,y E B{0,1 — e); 
(Hi) Furthermore, for even v 

(2.16) \Fi^y,-r,XT{t)'^{r^hD^i}u\ < C{T + l)h' Vr < e VT > \ogh\ 
and therefore 

(2.17) |</?(-e/iL)2)e(x,y,0)| < C/i^ 

Easy proofs are left to the reader. Inequalities ()2.14j) and ()2.16|) are proved by the 
standard elliptic methods. Then inequalities ()2.15|) and (|2.17|) are proved by the standard 
Tauberian methods. 

Thus in what follows we localized u with respect to hD^^ (and thus with respect to 
hDy^ as well due to symmetry. This localization is more precise than Xi,yi localization as 
e{h\ loghl)'" < ji < Ch~"; as yU < e{h\ \ogh\)~" both localizations are equivalent. 

2.2 Outer zone 

In this subsection we consider outer zone 2^out = {C*o7 < < 7i} with 

71 = min(e,C(/i/i)-^/("~^)). 
This definition works fine under ()2.8p but should be modified as 

(2.18) e{h\\ogh\)-^ < 11 < CqH-^ 
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to 

(and in addition {^i > Co} as u is even)^-'. 

Our first statement is that in the quantum dynamic in this zone the magnitude of xi 
(and ^2) persists. 

Proposition 2.6. (i) Let x E 5(0, |) fl {x,±Xi > C7}, 7 = e|xi| satisfy p. 6|) and E 
C^{B{x,'j)) be a rescaling of the standard function. Let ipi he ''y -admissible and supported 
in {l/(2Co) < xi/xi < 2Co} and equal 1 in {I/Cq < xi/xi < 2C07}. Let Ti = e/i7^. Then 
for T = Ti 

(2.19) iFt^h-wXrmi - Ax)uiJy\ < Ch^ 

(a) Let ^2 = ±Pj P > C and (p E C^{[j, |]) and (pi E C^, ipi = 1 on [|, |] he the standard 
functions. Let Ti = ep. Then for T = Ti 

(2.20) \Ft^h-wXT{i)U-MThD^))MThD,)\^C'^'- 

Proof. Note that both statements are equivalent as p x fi'y'^ so we will prove (ii). As > 7 
we constructed by p.27|) symbol which after rescaling becomes X = ^2 ~ P~^x\~'^ P{x)pi 
such that 

(2.21) \{a, X} - aa\ < Co 

Then using arguments of the proof of theorem 3.1 [Brlvrj with symbol 

(2.22) zu{<iCt±X) 

where w is the same function x as in this proof and ^ = ±1, one can easily prove (ii) as 
long as symbol ()2.22j) is quantizable; this condition is equivalent to ()2.6p . 

Surely one needs to check that in time Ti we stay in -8(0, 1) but it will be done in the 
next proposition. 

If (j2.6p is violated (and thus (j2.8p is violated as well) we can quantize with respect to 
(^2,4^2) but not with respect to {xi,^i) unless Xi > Ch\ logh\ and instead we note that 

{a, ^2} = aa + a^ [^2 - 4>Qu/ +0(1) 

and then X = ujS,2 with uj = (j){0, 0:2)"^ satisfies (I2.2ip and therefore the same arguments of 
the proof of theorem 3.1 |BrIvrj remain valid but we consider h-pdo with operator- valued 
symbols (in the axillary space K = L^(M^)). □ 

^^Restriction to |xi| will be |xi| < Ch\\ogh\. 
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The following proposition estimates by Cqu'^j'" and Cop~^ X2-speed of the propagation 
from above; the same result will hold for inner zone as well. 

Proposition 2.7. (i) In frames of proposition [KW i) let ip he supported in i-vicinity of y 
and ip2 = ip2{x2) satisfy 

(2.23) ^2 = 1 m {\y2-y2\<Co{^lYr'T + Co^] 

with i = Cofx-^-f^-" . 

Then for Tq = Coh\ log h\ < T < Ti 

(2.24) |F,^;,-i.XT(t)(l - i^2.)u^|Jy\ < Ch' Vr G (-e,e); 

(a) In frames of proposition lKW ii) let if) = '0(3^2) be supported in £-vicinity 0/I/2, (p = ^{^2) 
in ep-vicinity of ^2 and ip2 satisfy ^2.2^ . Let 

(2.25) lp>Ch\\ogh\. 
Then for Tq <T < Ti 

(2.26) \Ft^h-^^XT{t){^-^2x)u^y^\hDy,)\ < Ch' Vr G (-e,e). 

Proof. Both statements are equivalent under condition (j2.6j) . 
(i) Assume first that ()2.6p holds. Then symbol 

(2.27) Y = X2- 0- V'a;}-"^! 
satisfies 

(2.28) I {a, y} I < Cop-' 
with p = p'j'^ and symbols 

(2-29) w{Co^^±-^) 
{<; = ±1) are quant izable as long as 

(2.30) T > To = C(7/i|log/i|)^ 
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and then standard energy estimates method (see theorem 3.1 |BrIvr j ) imphes ()2.25j) with 
ip, and ■ip2 replaced by h-pdo supported and equal 1 respectively in domains 

{{y,v) ■ \yi - yi\ < \y{y,vi) - < 

: \yi - yi\ < €7, \Y{x,i,) - Y{y,i,)\ < ^ + C/i-^T}, 
as To < T < Ti = C/i7'^ and 

(2.31) 4 = C(7/i|log/i|)V"^ < ^ < e- 

Original estimate ()2.24|) follows from this because \Y{x,^i) — X2I < CqE. 
As fl2.24|) is proven for T = Tq it is valid also for T G [Tq, Tq]. 

(ii) Assume now that ()2.(ij) is violated. Then xi is localized only in Ch\ log vicinity of 
rather than in 7-vicinity of it and Xi is not separated from 0. However a naive idea that 
^4>Xi/h' should be replaced by ^2 and thus we should define 

(2.32) Y = X2- ly-'xiU^' 

instead of ()2.27|1 works; (j2.28j) and ()2.29j) are preserved, also (j2.3()j) becomes Tq = Ch\ log/i| 
and (jOT]) becomes (jT^ . □ 

The following proposition estimates by eo/i~^7~'' (and almost equivalently) by €qp~^ 
X2-speed of the propagation from below; this result would fail in the inner zone unless far 
from periodic trajectories. 

Proposition 2.8. (i) In frames of proposition \2. (k i) let iIj he supported in -vicinity of y 
and il)2 = ip2{x2) satisfy 

(2.33) ^^2 = 1 in {\y^-y,\>eo{firr'T-Coi} 

with i = fi'^-f^-". 

Then for Tq = Co7 < T < Ti inequality <\2.2^ holds. 

(ii) In frames of proposition Wl^ ii ) letip = ip{x2) be supported in i-vicinity 0/I/2, = '^{^2) 
in ep-vicinity of ^2 o-nd ip2 satisfy ^2. . Let i satisfy 1^2. 2(^ . 
Then for Tq = Coh\ logh\ < T < Ti inequality !\2.27^ holds. 

Proof. One needs to apply the standard arguments with 
(2-34) w{Co,^±^) 
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assuming that inequality ()2.29|1 is reversed to 

(2.35) |{o,F}| >eop-^ 

In frames of ()2.(ij) one needs to modify Y in the spirit of p.28j) which would be equivalent 
to the following modification in frames of (ii). Namely, note that as 7 ^ 1, p ^ 1 

{a,Y} = ^^'{-u-y,+pl) + o{p-') = i(l + u-')^2\pI - pD + ^(1 - '^~')^2\pI + pD 
and if we redefine Y as 

(2.36) Y = X2- v-^xiiii^^ + ]^v-\l + z/-^)xieiM2"' 
we arrive to 

(2.37) {a, Y} = \{1- u-')i^\pl + p\) + o{p-') = 1(1 - u-')^^\a + V) + o{p-') 

which implies fj2.36|) and also symbols ()2.35|) are quantizable. □ 

From the proof of proposition 12.81 immediatelv follows that ()2.2(ij) holds with iIj{x2), 
ip2{,X2) replaced by ip{Y{x,^))^ , ip2{X{x,^)Y with assumption £ = Cop^^j^"'^ replaced by 
flf^ i.e. 

i = Cop-^-f-''h\logh\. 

Then we immediately get 

Corollary 2.9. (i) In frames of proposition W^ i) 

(2.38) \Ft_,h~^rXT{t)T{u^) \ < Ch' Vr G (-e,e); 
as Tq = Ch\ log/i| <T <Ti and thus 

(2.39) Wt^hMxrAt) - Xfo(i))r(MV;)| < Ch' 
where I remind that x is supported in [—1, 1] and equal 1 in [— |, |]. 

(ii) In frames of proposition }^. bV ii) ^2.3^ . 1^2.39^ hold with ip{xi, X2) replaced by ip{x2)f{hD2) 
with (f supported in {\C,2\ > C}. 
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Assume now that ()2.6|) is fulfilled. Then we arrive to 

(2.40) iFt^H-wXTitm^n < C-fh~' VT G [To,Ti] 

as ip = ip'{xi)tp'(x2) is an element of 7-admissible partition with respect to xi] really this 
estimate holds with T = Tq due to the non-degenerating results rescaled ( |Ivr6j ) : x t— > x/7, 
h ^ h' = h/j, t I— > t/'-f and /i hh> ^' = fij'^ with fi'h' < 1. Then contribution of each 
element of 7-admissible (sub)partition with respect to x contributes C{h/'-f)~^ x 7 with 
extra-factor 7 coming from dt in Fourier transform Ft^^-ir- After summation with respect 
to all X 7~^ such subelements we get (j2.4(J|) . 

After ()2.4()|1 is established we can apply Tauberian arguments and conclude immediately 
that the contribution of ip to the remainder estimate does not exceed 

(2.41) CTf ^ X /I'S = C(/i7'')-i/i-S = C/i-^^"'' 
while the principal part is given by formula 

(2.42) h'^ f (Fi^,-i.XT(t)r(M^))rfr 

J —00 

with arbitrary T GG [To,Ti]. 

Then summation of (j2.4ip over zone |a;i| > 7 results in C^^^h"^'-/^^'^ coinciding with 
C'yh~^ as 7 = 7 and we arrive to 

Proposition 2.10. Under condition 1^2.8^ contribution of the whole outer zone Zont = {x : 
l^^il > C'07} to the remainder estimate 

(2.43) n = \rie^) - h-' f (Ft_,^,-.,XTit)T{u^P)) rfr| 

J —00 

does not exceed C'yh~^ where ip = 1 — (j){xi / {Co'j))ip' {xi)ijj" {X2) with all functions regular, 
uj supported [—1, 1] and equal 1 on [— ^, 

The similar results as condition ()2.8|1 violated will be derived later in sections 4, 5. I 
will replace ()2.42j) by more explicit expression in section 3. Basically it will be the answer 
prescribed by the non- vanishing magnetic field theory but for /i large enough some correction 
will be needed. 
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2.3 Inner zone. I 

Inner zone is Zjnn = {\^2\ < Co} or under condition ()2.8|) equivalently (around energy level 
0) Zinn = < Co7}; under condition ()2.18|) |xi| < Coh] log/i| in the microlocal sense in 
this zone. For a sake of simplicity of notations we use 7 instead of 7 until the end of the 
section. 

Let us start from the confinement of propagation in the inner zone and a finite speed 
of propagation; proofs repeating those of the corresponding propositions of the previous 
subsection are left to the reader. Just to remark that in this zone we study ^2 and X2 
directly without going to X and Y. 

Proposition 2.11. fi) Let condition p.<S|) be fulfilled. Let x G B(0,^)r]{x,\xi\ < 7} and 
ip G C^{B{x,'~f)) be a resettling of the stttndttrd function. Let ipi be '-/-ttdmissible, supported 
in {\xi\ < 3C07} ttnd equttl 1 m < 2Co7}. Let Ti = e. Then estimttte ^2.19^ holds 
with T = Ti . 

(11) Without condition {13) let^ € C^{{\, \)) , if £ {{-Co,_Co)) , ifi e {{-3Co,3Co)) 
ttnd equttl 1 in (— 2Co, 2Co). Then estimttte 1^2. 2(J^ holds with ^2 = 1 o,nd T = Ti. 

Proposition 2.12. (i) In frttmes of proposition lKTW i) let = ^{x2), 

4^1 = 1 in {\y2 - ml < CqT}. 
Then for Tq = C07 <T<Ti inequttlity 

(2.44) |F,^;,-i.XT(t)(l - ^Pi.)uQy\ < Ch' Vr G (-e,e); 
holds with Qy = ipy. 

(a) Without condition (j^./yjl this estimttte holds with Qy = il){y2)<^{—hDy^) with ip G 
Co°°((-Co,Co)). 

Now we are studying the finite speed with respect to .^2 and the proof repeats those of 
proposition 12.71 

Proposition 2.13. Let^) = ^{x2) G C^{{-\,\)), ^ G C^{B{f]2.p)) ttnd^i G Co°°(B(r/2, 3p)) 
is equttl to 1 in B{f]2, 3p) be stttndttrd functions resettled. Let 

(2.45) p>Cih\\ogh\+Ca 

^^Cf proposition 12. 61 
6)Cf proposition O 
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Then for Tq = C27 < T < Ti (p) = ep inequality 
(2.46) \Ft_^h-wXT{t){l - ViihD,,))u4jy^{-hDy,)\ < Ch' Vr G (-e,e 



We will make more specific statements later. Now we want to prove that as ^2 is disjoint 
from k*V'2 there is a drift with the velocity x (^2 — k*V'2) (with the correct sign). Let 
V = ^{^2} £ C^(^B{f]2, p)) and = 4'{x2) G C^{B(jj2,tj, i = eip be standard functions 
rescaled. Assuming that 

(2.47) p>C{h\\ogh\)-2 
we conclude from propositions 12.121 12.131 that 

|F,^,-i.XT(t)(l - MhD.2)Mx2))utlj{y2M-hDy,)\ < Ch' Vr G (-e,e) 
as long as T < ep. 

Proposition 2.14. Let conditions and fTOTI he fulfilled. Letip = 'ip{x2) G C^[B{y2, i)) 
and (p G C^(-B(f/2, p)) also standard function rescaled, with p < Cq, 

(2.48) T]2-k*V{0,y2y^ = ±p. 

Let tpi = ipi{x2) be also a standard function rescaled and equal 1 in {\x2 — y2\ > ^opT — Coi} . 
Then estimate i\^-44\ l holds with Qy = ip{y2)'p{—hDy^) as |r| < ep, i > 7, 

(2.49) pi>Ch\\ogh\, 

(2.50) T > Co-. 

P 

Proof. With no loss of the generality we can assume that V{0, 1/2) = 1. Let us consider 

(2.51) Y = X2-7{-,^i) 

7 

where Z is a function introduced in proposition 11.61 and pick up 7 = exactly. Then 

I {a, Y} — P{^2 — k*) — aa\ < e'p with a constant e' which one can make arbitrarily small. 
Then we can apply standard arguments with 

(2.52) w(Co<;]{Y{x,0 - Y{y,v))) ± ^) 

as long as Tp > Ci; symbol (j2.52p is obviously quantizable. □ 
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It follows from the proof that estimate ()2.44|1 holds with Qy = tp{Y)'"(j){—hD2) and tpi 
replaced by ipilY)"" without assumption i > '~f] this estimate immediately implies 

Corollary 2.15. In frames of proposition estimate ()^.56|1 (see below) holds with the 
standard function ip and |r| < ep as long as 



(2.53) P>Co{h\\ogh\)-2 

i\] 
P 



(2.54) To = <T <Ti = ep. 



Note that Tq < e7 as long as 

(2.55) p>p/=^C(-)^ 

7 

and we arrive to 

Proposition 2.16. Let conditions (fF^ . (fTTJj) and he fulfilled. 

Then 

(i) Estimates 

(2.56) |Fi^;,-i.XT(t)r(<^(/iD2)H I < Ch\ 
and 

(2.57) |Fi^,-v(xT,(t) - Xfo(i))r(¥^(/iA)V^u)| < Ch' 

hold with \t\ < ep and Tq <T <Ti = ep. 
(a) Furthermore 

(2.58) \Ft^h~wXTAt)^ivihD2)iPu)\ < Ch^'^p. 

Proof. We need to cover T < e'j only. After rescaling x ^ x/j, t ^ t/j, h ^-^ h = h/'j we 
find ourselves in frames of the standard propagation and therefore estimate ()2.38|) with h 
replaced by h would hold as long as Coh\ \ogh\ < T/7 < e. Then if h < or equivalently 
/i < h^~'^ with arbitrarily small exponent 5 > we arrive to original (j2.38|) with indicated 
T, r. 

Further, these arguments imply ()2.58p with no such restriction. 

For h^'" < p < e{h\logh\)~'^ we instead can make a fine (e'7, e')-subpartition with 
respect to (a;i,^i) and then depending on subelement apply standard arguments with one 
of 

(2.59) w(Co^^(xi - y,) ± |), w{Co^^i^^ - Vi) ± ^] 
which are quantizable. Easy details are left to the reader. □ 
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Then applying Tauberian arguments we conclude that contribution of this element to the 
remainder estimate does not exceed Ch"^^p/Ti x Ch~^^. After summation over partition 
'm p > ei we get remainder estimate Ch~^'-f] we will cover smaller values of p in the next 
subsection. Thus we arrive to 

Proposition 2.17. Under conditions contribution of zone 

{kil < 7, 16-^*^(3:2)^1 >eo} 
to the remainder estimate does not exceed C^p^^^'^h^^ . 



2.4 Inner zone. II 

Now we need to consider "near-periodic" zone 

{\xi\<^,\^2-k*V{x2)--\<eo}. 

Let us investigate this case under condition ()1.2H) . assuming first that d^^V is disjoint 
from 0. We consider a bit more general case: namely dx2V x ( with large enough parameter 
(. We consider first propagation assuming that that l^ajj^l < (. 

Proposition 2.18. Let conditions 1^2.8^ . be fulfilled and 

(2.60) Id^^V] < C as xi = 
with 

(2.61) C > Co7- 

Let ip = 'ip{x2) G C^{B{0, |)) be a standard function and ip = </?(6) be a standard function 
rescaled supported in {\^2 — k*V^{0)\ < p} '^^ with 

(2.62) p > CoC- 

Let ipi = v^i(6) be also a standard function rescaled, equal 1 as |^2 — k*V^ {0, X2)\ < Sp. 
Then 

(i) Estimate p.y/6D holds as |r| < eC, and ChC,~^ \ \ogh\ <T<e; 
(a) In addition, if 

(2.63) p>pi = iC^-^h\\ogh\)-^ +C^. 
■^^Then under condition IC2 - k*V^O, X2)\ < 2p for all X2. 
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estimate (\2.4(^ holds as \t\ < £( and 

-1 

(2.64) C hC^ \ log h\ <T <€mm{-,^). 

Proof. I leave to the reader a standard proof based on the axillary symbol 

(which under our assumptions is quantizable) that the speed of propagation with respect 
to ,^2 does not exceed C(. This implies (i). 

To prove (ii) one needs to show in addition that the "averaged" speed of propagation with 
respect to X2 does not exceed p; under condition ()2.64|) one can easily do it for t = T*V{OY^^ 
by means of the standard microlocal analysis. I leave details to the reader. More delicate 
analysis of subsections 2.6-2.9 will also imply results of this subsection. □ 

It immediately implies 

Corollary 2.19. Under conditions 1^2.1^ . 1^1.21^ . ^2. 6U^ — ^2. propagation remains con- 
fined to 

{\x2\<\M2-k*V{Q)^<p + C + m 

as T < e min(p~^, C~^) ■ 

Now let us estimate from below the propagation speed with respect to ^2- 
Proposition 2.20. Let conditions (fF^ . (f07|) . (fE^ he fulfilled and 

(2.65) eoC < \dx^^\ < C as Xi = 

with ( > Co7- Let ip = ^{^2) be a standard function rescaled supported in p-vicinity of 
f]2 = k*V{0)2. Then estimate !\2.4(^ holds as |r| < eC, 

(2.66) To = CpC' <T<T, emin{p-\ 
= 1 - V?. 

Proof. I leave to the reader a standard proof based on axillary symbol 

^{Co<;^{^2-V2) ± ^) 
which under our assumptions is quantizable. □ 
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Then in frames of this proposition ()2.56j) holds. Note that Tq < e'j as 
(2.67) C> C-f^^h\ log h\ 

and due to the same arguments as in proposition 12. 161 we arrive to 

Proposition 2.21. Let conditions of proposition W^W\ and \2. 61\ he fulfilled. Then esti- 



mates dOTI), (fF:^ hold as To = Ch\ \ogh\ <T <Ti with Ti defined by (fF^ . 



Then due to Tauberian arguments we arrive immediately to 

Proposition 2.22. Under conditions ()i^.<^ . (ji.i^ij) and 6'7D the contribution of the inner 
zone to the remainder estimate IZ does not exceed C n^^^'^h"^ . 

In particular C ^ 1 satisfies (j2.67|) under condition (j2.8j) . Case of smaller C, < Ch'~f~^ \ logh\ 
will be considered in subsection 2.7. 

2.5 Inner zone. Ill 

Let us finish analysis in the "near-periodic" but not "periodic" zone. To do this I am 
going to analyze elements of (0:2, ^2) partitions on which either |^2 — V^^'^k*] + \dx2V\ > Cpi 
branding the rest as "periodic zone" . 

The following results would be the results of the previous subsection rescaled but some 
conditions are more relaxed because after rescaling xi— 71— >7/£. C ^ our new 
symbols have smaller derivatives (with a factor €) with respect to xi than it was assumed 
there. 

Proposition 2.23. Let conditions ^2.^ . ^1.21\ . ^2.6S^ be fulfilled. Consider point {y2,V2) 
and assume that 

(2.68) \r]2-k*vH0,y2)\<p, 

(2.69) \d;^^V\ < C„C^^-I"I as xi = 0, \x2 - ml < i Va : |a| < K, 

(2.70) ip>Ch\\ogh\, 
1 

2' 



(2.71) ec < Ip, 



(2.72) i>Coi, P>Co^, C>Coi + Coh 



*)Cf. propositioiEU and corollary inn 
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Then for \r\ < p, 

(2.73) |t|<r/=^^min(^,^) 

quantum evolution starting from {1^2 — ^2! ^ ^^AV2 ^ ^ |p} is confined to {\y2 — y2\ ^ 
Proposition 2.24. Let conditions of proposition {2~2M be fulfilled. 

(2.74) \r]2-k*vH0,y2)\<P, 
then estimate ()i^.56|) holds for \t\ < p, 

(2.75) To''='C^^^^<T<n; 

P 

(a) If (\2. 6^ is fulfilled as xi = 0, |x2 — ^2! < ^ t/ien estimate ^2.56} holds for |r| < p, 

(2.76) To"=^C^^i^<T<Ti. 
Now, given point (1/2, r/2) let us define 

(2.77) ^ = p = C = Qim. m), Qix2, 6) = e (I6 - ^ V(0, ^2)2 1 + |a,,^(0, X2) |) . 

Note that (jTZU)), (jT7^ become equivalent to (PIk?|) and thus propositions 12221 hold 
with Ti = e under condition ()2.63|) . Also note that Tq < 67 is equivalent to condition ()2.63|) . 

Then using the same arguments as above we get ()2.56j) with To = Ch \ logh\ < T < Ti 
and also estimate 

(2.78) \Ft^h-wXTit)^HhD2)iJu)\ < Ch~'pji Vr : |r| < ep 
and Tq <T <Ti. 

Therefore due to Tauberian arguments contribution of (£, p)-vicinity of (2/2, r/2) to the 
remainder estimate does not exceed the right-hand expression of ()2.78|) multiplied by CTf ^, 
i.e. Ch~^'~fpi] I remind that i = p here. Then for given p the contribution of zone 
{{x2,(,2), \p < Q{x2,i2) < p} to the remainder estimate does not exceed Ch'^pj. Finally, 
summation over g results in Ch~^'j. This yields statement (i) of the following proposition: 

^)Cf. proposition 12201 
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Proposition 2.25. Under conditions ()g.<^ . 1^1.21^ 

(i) Contribution of zone {{x2iS,2) '■ q{x2,C,2) > ^i} to the remainder estimate ( 2.4^ with 
any T E [To,e] does not exceed C fi~^^'^h~^ ; 

(a) The total remainder estimate l \2.43^ does not exceed 

(2.79) C/i-^/i-i + Cpi/i-i; 

(Hi) In particular, as pi < 7 which is equivalent to 

(2.80) fi < C{h\ log h\y^ 

the total remainder estimate does not exceed C^~^^'^h~^ . 



Proof. To prove (ii) one needs to estimate the contribution of periodic zone {g < Cpi}. 
One can take Ti = e7 tliere and tlius its contribution does not exceed C/i~^pi7 x 7"^. □ 

2.6 Periodic orbits. I. Pilot-model 

Therefore one needs to consider contribution of the periodic zone 

(2.81) Zper = {(X2,6),f?(a:2,6) < ^1 = C{^''h\hgh\)"^} 
as ()2.80|) is violated i.e. as 

(2.82) C{h\logh\y^ <fi< e(/i|log/i|)-^ 
Assume that proposition 12.241 remains valid^^-* as 

(2.83) Po = C{h\ log h\y^ + C-f<g<pi. 
Then for particular partition element inequality 

(2.84) \Ft^h-irXT{t)r{^ihD2Mx2)u)\ < Ch~'g^^ x ^^^^ = logh\ 

holds as |r| < ep, Chp~^ logh\ < T < Ti = e where the second factor in the middle 
expression of ()2.84p is just Tq. Inequality ()2.84p leads to the estimate of the contribution 
of the given element to the remainder, equal to the right-hand expression of ()2.84jl (since 
Ti X 1). Then the total contribution of zone {(^2,^2) '■ \p < q{x2^^2) ^ p} is equal to 

Actually it does but with a twist. 
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Ch~^'-f p~^\ log h\. This estimate is not only is not good enough for the sharp remainder 
estimate but is much worse than estimates we had before^^-*. 

To improve the estimates above let us start from the pilot-model operator. However let 
us consider first 1-dimensional operator 

(2.85) L{k) = ^ {h^Dl +{k- x\lvf - l) 
with Yi= h/'^. 

Proposition 2.26. As h<^l 

(2.86) Spec(L(fc)) n [-e, e] = ^(^(^ + \)k k, h) + 0{h^), neZ 

and consists of simple eigenvalues where g{T] k, h) is an analytic function of all its arguments 
as \t\ < e, h < e,\k — k*\ < e, inverse (with respect to r) to 

(2.87) g{T; k, h) = T{k)~^ j T{k{l + 2r)~5) (i + 2r)^ dT-hi + 0{h^) 

where T{k) is an elementary period of section 1 calculated as = 1, r = 0. 
In particular 

(2.88) gift-iT(%(L;fc,a) ^ Xe{L)I. 

Proof. Note that for L{k) Hamiltonian flow on the energy level r is periodic with period 
T(r, k) =T(fc(l + 2r)-^)(l + 2r)^. 

Proof now follows from Bohr- Sommerf eld formula. To derive asymptotics with 0{fii^) 
error one should note that subprincipal symbol of L{k) is and Maslov' index lm of the 
trajectory in (a;i,.^i) space is 2. 

□ 

Proposition 2.27. Let condition i\2. 8^ be fulfilled and let Q = ip{hD2) he a partition 
element supported in B{k*,p). Let us assume for simplicity that V{0) = 1. Let A be the 
pilot-model operator: 

(2.89) A = A = ^ (h^Dl + {hD2 - px^/uf - l) 

One can increase Ti and decrease po in the pilot model cases improving this estimate but H2.84|l should 
be improved as well. 
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(a) First pilot-model, v ~2. k ~ 0.65 (b) Second pilot-model, v ~ 2^ 

k = 



Assume that 7 > Ch\ \ogh\. Then modulo negligible operator 

(2.90) e''-'^'*-^QUA) ^ e^'-'^'^QUA), 

(2.91) A = g{A,k*,h), T* = T{k*), 

B = /3{hD2, A, h) such that 



(2.92) /3(6, r, h) = «:(6, h) - K)' + K,{h) + 0(|r|) 

and function /3(^2,t, 0) coincides with b defined by l \1.4^ - Further, = k* + 0{h), Ki = 
0{h) are analytic functions as h 1. 

Proof. Proof immediately follows from propositions 11.151 12.261 □ 
Remark 2.28. From now in the analysis of periodic zone k* = k^ and respectively p = 

\^,-kl\ OT p = - klV^'l 

Now (|mn|) implies that 

(2.93) e''-''-^Qx.M) = e''-"'"^e''"''''QXep{A), t' = ^[-\,t" = t~ t' 

7 

as \t\ <Ti = ep^^. 
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Here important thing is that on the partition element in question operator e^'^ is a 
proper FIO as long as 

(2.94) |t'|<T{ = -^ 

which is greater than Tq = Cp~'^h\ logh\. Then repeating corresponding analysis of subsec- 
tions 2.4, 2.5 one can prove easily 

Proposition 2.29. Let p be defined as in remark V2. Then proposition \2.1t^ ii). corollary 
\2.iy\ and proposition]^. 24\ remains valid with ^2. replaced by <^2.8!^ . 



These statements extended allow to prove 



Proposition 2.30. In frames of proposition 2.21 assume that supp(v9) fl B(k^,p/2) = 
and that ip is a fixed admissible function. Then 

(2.95) \Ft_,h-r,XT{t)^{^ihD,)ij{x2)u)\ < Cp~'\\ogh\\ 

Proof. Due to proposition 12 .241 extended by proposition 12.291 it is sufficient to consider T = 
To only. Consider left-hand expression of ()2.95|) with Xr(t) replaced by Xn(t) = X-y(t~''T'T*'y), 
IM < Tq/T*'-j. Note that then e^^ "'^^ is /i-pdo with symbol which will be regular after 
rescaling xi ^ Xi/ (7I log/i|), ^1 ^ ^i/| log/i|, ^2 ^ G/| ^ogh\. Therefore 

(2.96) \Ft^n-.,Xn{t)T{v{hD2)i^{x2)u)\ = |F,_^-i.Xo(t)r (Qn^^) | 

with Qn = e"^ ^^^^ (p{hD2)ip{x2) and therefore expression (j2.96p does not exceed Cp'yh~^ x 
h\ log/i| = Cp'jh~^\ \ogh\ due to the standard theory. 

After summation by n : |n| < CT0/7 we get Cp7/i~^| logh\ x hp~'^'j~'^\ logh\ which is 
exactly the right-hand expression of ()2.95p . □ 

I believe that factor | log/ip is superficial and one can get rid of it: one factor | log/i| is 
due to an inaccuracy in the estimate of ()2.95|) modified and another due to the lost factor 
(1 + n7/f^)"' in the same estimate (with Tq = Cp "^h). However even such weakened 
estimate is sufficient for our needs unless p is really close to h^^ and the latter special case 
will be considered separately. 

Really, after ()2.95j) is proven, Tauberian theorem implies that the contribution of this 
partition element to the remainder estimate does not exceed Cp~^ \ log /ip x T]~^ = C\ log /ip 
and thus the total contribution of zone {|^2 — k* — fi\ > C7} does not exceed this expression 
integrated over p~^ dp resulting in C\ log/ip while contribution of zone {|.^2 — k* — h\ < C7} 
does not exceed C'-fh~^ as we already know. Thus we arrive to 

Proposition 2.31. For operator coinciding with the pilot-model i\2. 89^ in B{0, 1) remainder 
estimate^"^^ is Cjh~^ as p < Ch~''\\ogh\~'^'' . 

i2)which still means TZ given by (|2.43(l . 
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2.7 Periodic orbits. II. Another example 

Consider now the pilot-model ()2.89|) perturbed by potential V{x2) satisfying 
(2.97) eoC < \d,,V\ < ( Vx. 

Then we can apply the same arguments as above for partition elements with > C thus 
covering completely case ( < 6*7^. On the other hand, as C > 7^ we replace p by ^ = 
max((^2^p) and on partition elements with g = we again can apply the same arguments 
with Ti = eC~2, Tq = eh(~^ \ logh\. 
So, we arrive to 

Proposition 2.32. For operator coinciding in B{0, 1) with pilot-model perturbed by po- 
tential V{x2) satisfying <\2. 9(J^ remainder estimate does not exceed 0{^h~^) as long as 
+ 7 > Ch\\ogh\^ i.e. either p < C/i"^! log/ij-^^ or C > Ch^\ log/i|^ 



2.8 Periodic orbits. III. General case 

Consider now general case. Let V9(^2) and ip{x2) are two p- and ^-admissible functions 
supported in p/3- and £/3- vicinities of ^2 and X2 respectively. Let us assume that conditions 
(frH8|l - (j2T2|l are fulfilled. 

Then due to proposition 12.241 extended by proposition 12.291 as 

(2.98) |t| <T/=i:'emin(^,^) xep£(p2 + £C)"' 

propagation started from suppy? x supp?/; is confined to (p/3, £/3)-vicinity of it. 
Now we can easily generalize proposition 12.271 

Proposition 2.33. Let conditions (f^ . (fTT?|) . - ( [E7^ be fulfilled. Assume for 

simplicity that V{Q) = 1 

(i) Identity (f^^ holds as A is defined by jUSII); 
(a) Therefore (f03|) holds. 

Proof. Proof repeats one of proposition 12.271 I leave details to the reader. □ 

Proposition 2.34. Let conditions (IF^ . (fTT?|l . (fO^ - ( [E7^ be fulfilled. Assume for 
simplicity that V{0) = 1. Then estimate 

(2.99) \Ft^h-wMt)'^{'fihD2Mx2)u)\ < CTpih~'\ \ogh\ 
holds as 

(2.100) 7 < T < = Cih\ log/i|(p2 + 
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Proof. Let us note that for \t'\ < Tq operator e*'' remains legitimate /i-PDO and there- 
fore contribution of the time interval [^1,^2] with \tj\ < Tq, \t2 — ti\ = e7 does not exceed 
Cpi'jh~^ \ logh\ and therefore the contribution of x T/7 of such intervals does not exceed 
the right-hand expression of ()2.99p . □ 

It immediately implies 
Corollary 2.35. Let in frames of the previous propositions either 

(2.101) \^2-klV{0,X2)\>p, 
or 

(2.102) \d,,V\ > eoC 
in i-vicinity of X2 ■ Then 

(i) Left-hand expression 1^2.99^ with T <Ti does not exceed CTQpih~^\ log Zip; 

(a) Contribution of the partition element (f{hD2)ip{x2) to the remainder estimate does not 

exceed C\ log/ip. 

Proof. Statement (i) immediately follows from propositions 12.241 extended and 12.341 state- 
ment (ii) follows from the Tauberian arguments: contribution of this element to the remain- 
der estimate does not exceed C/i~^p^ | log /i| xTq/Ti = C| log /ip since Tq/Ti = Ch\\ogh\/{pi) 

□ 

Assume now that 

(2.103) \d,,V\ + \dlV\ > eo. 

Let us introduce functions g by ()2.76|1 and let us consider ^-admissible partition in (^2,^2) 
and apply the same arguments with p = C = i = Q] then on each partition element either 
\^2M — k'^V{0,X2)\ > eoQ or \dx2V{0,X2)\ > eg and everything works as long as 

(2.104) ^>C7, g>C{h\logh\)^ 

where the second inequality is equivalent to dp > Ch\ \ogh\. 

We know that the contribution of each such element to the remainder estimate does not 
exceed C\ log/ip and therefore the total contribution of all such elements does not exceed 
Clog/ipJ with / = / g-'^dx2d^2 ^ C\ \ogh\ due to condition (|2.1()3jl . 

On the other hand, as g < g = C7 + C{h\\ogh\)2 let us redefine g as g. Then 
contribution of each such partition element to the remainder estimate does not exceed 
Cg^h~^ = C'-f'^h"^ + C\ log/ip and again due to condition ()2.103|) there is no more than C 
of such elements. Thus we arrive to 
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Proposition 2.36. For operator satisfying in -8(0, 1) conditions 1^1.21^ and 1^2.10,'^ the 

remainder estimate is 0{'-fh~^) as long as fi < Ch~'^\ logh\~^'^ . 

On the other hand, exactly the same approach but without condition (|2.103p results in 
the remainder estimate 

C g-^\ log h\^ + Cgh-^ 

where the first term is an upper estimate of C\ log/ipJ with integral over {g > g} and the 
second term is the contribution of elements with g^ g. Picking up ^ = C max(7, /la | \ogh\) 
we arrive to estimate 

C-ih^^ + Ch~^\ogh\ 

which is our target C'jh~^ as long as 7 > /i^j \ogh\ i.e. /i < Ch~'^\ log/il^*^. Thus we arrive 
to 

Proposition 2.37. For operator satisfying in B{0, 1) condition 1^1.21^ 

(i) The remainder estimate is 0{'yh~^) as long as fi < h~'^\ loghl^"; 

(a) The remainder estimate is 0{h~^ logh\) as long as h''^] loghl^^fi < Ch^^\ log/ij^'". 

2.9 Periodic orbits. IV. General case (continuation) 

Finally I am going to prove the most general 
Proposition 2.38. Under conditions p.<j^ . () and 

(2.105U lC^I>^o. 

l<fc<m, 

with W{x2) = y(0,X2) the remainder does not exceed C{\\oghY^^^ +7/1""^) and thus the 
remainder estimate 0{{jh'~^^ as long as fx < Ch~'^\ logh\~^"^^^^'^ ; 

(a) Without condition \(2. 105)~^ the remainder does not exceed C (7"^^ +7/1"^) and thus the 
remainder estimate is 0((7/i~^) as long as fi < Ch~'^~^^^ with arbitrarily small Si > 0. 

Proof. Let us consider the scaling function 

(2.106U ^(X2) = 6i( Yl \d',W{x,)\^)^ +e 

l<fc<m-l 

with small enough constant ei > such that l^^^j^l < ^. Let us consider ^-admissible 
covering of (—1,1). 
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Then 

(2.107) \d^^,W\ < Coe{x2r-'' VA; : A; = l,...,m- 1; 

further, on each element with i > 2i one of these inequalities could be reversed with Co 
replaced by 1: 

(2.108) fc \d',^W\<Coiix2r-'' k = l,...,m-l. 

Furthermore, under condition |(2.105)n^ there is no more than Co of ^-elements with 
i G (/, 21) for each / and thus no more than C| log£| elements in total while without condition 
|(2.105)m] the corresponding estimates are C/~^ and Ci~^ respectively. 

On each ^-element oscillation of W{x2) does not exceed C£"^ and thus we can ponder 
condition 

\^2-KWiyy^\^p 

with p > Ci"^ and y the center of this element. As this condition is fulfilled and also 
pi > Ch\ logh\, p > h^~^ we can apply corollarv 12. 351^ 11) and conclude that the contribution 
of this element to the remainder estimate does not exceed C| log/ip. 

Now we need to sum contributions of all such elements. To do it in the most efficient 
way we redefine 

(2.109) ^ i{x2,^2) = ei{\^2-KW'^{x2)\+ KWi^2)\^)^ 

l<k<m-l 

Then for each given p the total contribution of all elements of this type would not exceed 
C|log/ip and Cp~™ | log /ip under condition |(2.105)^ and without it respectively; after 
integration with respect to p~^dp we get C| log/ij"^ and Cl~^ \ log/ip respectively and both 
these expressions do not exceed stated in (i),(ii) remainder estimates^^-* as we pick up 

(2.110) 1= Ci7^/(™-^) +Ci(/i|log/i|)^/("^+^). 

On the other hand, the contribution of each element with £ x £ does not exceed 
(^^m+i^-i I /;^|2 g^j^^ ^Yi.e total contribution of such elements does not exceed Ct^^^ \ log /ip/i" 



and Cf™"! log/ip/i ^ under condition (2.105)^ and without it respectively and these expres- 



sions do not exceed remainder estimates announced in (i),(ii) respectively. 

Thus we need to consider only elements with p < Ci"^ and i > 2£ with £ given by 
|(2.106)m] rather than |(2.109)m| Let us consider one such element {y — ^£,y + |£); after 
rescaling 

x^{x- y)/£, W{X2) - W{y) ^ {W{x2) - W{y))/£"', p ^ pr'" 



^■^^ Provided m is large enough in (ii) which we assume. 
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we find ourselves in frames of (2.105)m,_i- 

So in the rescaled coordinates let us define scaling function i[ by (2.109)m_i and redefine 
W as i'^^W. However, as we return to the original scale this new rescaling function 
ii = £[£ will be defined by (2.109)m-i with original W. Now we can repeat the same 
arguments as before with ii defined by ()2.110|) with m replaced by [m — 1). However now 
in the total partition there is no more than Ci of ii elements with ii G (/, 21) if original 
condition |(2.105)^| was fulfilled and no more than Cil ^ of them otherwise because as I 
mentioned after rescaling of the elements in question (2.105)^-1 always is fulfilled. Then 
we eliminate elements of ii partition with £i < 2ii and with p > Ci£^~^ estimating their 
total contribution by C|log/ip+^ + C-fh-^ and C\\ogh\"'+^ + Cl'^^H^^h"^ in frames of 
(i), (ii) respectively which do not exceed remainder estimates announced there. 

We will continue this process as long as our scaling function is defined just by 
(2.109)m-fc i-e. until k = m — 2. Then all elements remaining are eliminated by corol- 
lary under condition |9^.2W^| x C- D 

Thus I have covered the case 
(2.111) fi< h^^l log h\~^^ 

leaving the case of superstrong magne tic field /i-'^l log/i|~^^ < /i < C/i-^ for section 4. 



3 Calculations 

So far we recovered good remainder estimates but the asymptotical expression (|2.42j) I 
derived was implicit. In this section I am going to rewrite it in more explicit form as either 



(3.1) j S^^{x,0)ij{x)dx 

which is of mag nitude/i-2(l+/i/i) ^/^'^ or the same expression but with a smaller correc- 
tion term. Again I consider only the case (j2.111|) leaving the case of superstrong magnetic 
field for section 4. 

Remark 3.1. With no loss of the generality one can assume that = ip2{x2)'ipi{xi) with 
fixed admissible functions ipj; moreover, one can take ip=l as /i > Ch~^. Really, it follows 
from the analysis of subsection 0.1 that for fixed ip vanishing as xi = the remainder 
estimate is 0{'^h~^) while the principal part is given by p.ip . 
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3.1 Outer Zone 

I remind that the outer zone Z^ut = {C^ < \xi\ < 7i} and in this subsection 7 ranges 
from 7 to 7i as well. I also I remind that according to subsection 2.2 we can take here 
T = Ch\\ogh\ in flTT]) . 



Proposition 3.2. Let condition 1^2.8^ be fulfilled. Then 

(3.2) \h~' [ (Ft^f,-irXT{t)T{i:{l-i:,)u))dr- 

^^^(x, 0) {^{1 - ^1)) dx\ < Cifih + 1) 



for T = Ch\ logh\, fixed function ip and ipi{xi) = ipi{xi/'y) with a fixed admissible function 
ip, supported in (—2,2) and equal 1 in (—1, 1). 

In particular we get C'yh~^ -estimate as fi < Ch~'^ li'^^-^) , 



Proof. It follows from analysis of |Ivr6j that for 2-dimensional Schrodinger operator with 
non-degenerate magnetic field, binding parameter /j, and semiclassical parameter h under 
non-degeneracy condition |VV^/F| > eo the following asymptotic expansion holds 



(3.3) h ^ / (Ft^h-^^XTit)T{ilJu)jdTr^ ^ 

J —00 . 



-2+2n+2m 2m 



m,n>0 



as /i < /i^ ^ where 

(3.4) t^moh-^fih)^"" ~ J ^^^(x, 0)^{x) dx 



m>0 



and therefore one can replace ()3.3|) by ()3.4p with an error K02 + 0{{fihY). 

It also follows that in the same settings but with h^~^ < /i < Ch~^ another asymptotic 
expansion holds 



"^""^ m,n>0 m,n>0 



and comparison of this two expression with ()3.3|) where they overlap implies that k'^'^ = 
and an error in question is given again by K02 + 0((/i/i)^). 

Therefore in our setting the contribution (to this error) of any element of 7-partition in 
the outer zone will be of magnitude 1. To calculate the total contribution of the outer zone 
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one must integrate the latter expression over 7 ^^7 resulting in 0(7 ^) which is not as 
good as 0(7f^) with 7i = min(l, (/i/i) ^^^'^ which I claimed. However plugging /i7'^, 
h h/'j into 0{{jjihY) and integrating produces 0(7]"^) which is exactly what I claimed 
and one needs just to calculate explicitly the correction appearing from fi;o2-term. 

In the non-degenerate settings this latter correction does not depend on ix and therefore 
one can calculate it as /i = 0. But this is non-magnetic case and then due to the standard 
theory and condition \V\ > eo 



n>0 



(without integration with respect to dx) with xo(x)/i ^ = £^^(x, 0) and therefore the total 
contribution of this correction term is 0(1). □ 

Now I am going to derive a bit sharper remainder estimate; to do this I must appeal to 
the spectral projector of the corresponding one-dimensional Schrodinger operator: 

Proposition 3.3. Let condition (\1.21\ he fulfilled and let ipi be an even function, = 
ip{x2)- Then as h"" /(^I'-i) < ji < e{h\ \ogh\)~'^ , T = e'j the following estimate holds: 



(3.5) \h-' (Ft^h-wXT{t)^{^{^-i^i>))dT- 

J —00 

j (^^W(x, 0) - Sf'^ix, 0)) V^(X2)(1 - ^i(xi)) dx- 

j eo(xi,Xi;x2,6,0,/i)^(x2)(l - V^i(xi)) rfxrf^al < 

Ch~^^ + C\\oghf 

where eo(xi, X2, ^2, t, h) is the Schwartz kernel of the spectral projector eo(x2, ^2, t, K) for 
1- dimensional Schrodinger operator 

(3.6) ao(x2,e2,/i) = \{f^^Dl + {i2-x'i/uf - W^x^) 
with h = h/'j, while 

(3.7) S^'^ix, 0) = ^ V ^f2r + V{X2) - {2n + 1) fikx^'^) fih"' x^-' 

n>0 

is Magnetic Weyl approximation for corresponding pilot-model operator. 
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Proof, (i) To calculate u{x,y,t) let us apply first the method of successive approximations 
with unperturbed operator A obtained from A by freezing X2 = y2 V , a, (p: 

(3.8) A = i (^h'Dl + a(xi)2(6 - /i0(xi)xr)/z/)' - V^(xi)) , 

a = a{xi, ys), = 0(a;i, 2/2), = 'i^(a;i, 2/2)- 
Then the first term of 

(3.9) f (F,^,-i,XT(t)r(^/'(l - ^i)n)) dr 

J —00 

is also given by formula ()3.9|) but for operator A instead of A and it is 0(/i~^7iTo) = 
0{h-^^i\\ogh\) 14). 

The estimate for the second term gains the factor CT^/h x C/i| log/ip provided = 1 
identically or factor 

CT^{1 + /i7r'/i)//i X + 71^) I log/ip 

in the general case (I remind that then = 1 as Xi = 0); therefore the second term does 
not exceed Ch~'^^i {h + 7^) | log . 

Finally, in two-term approximation the remainder gains one more such factor and thus 
does not exceed C/i~^7i(/i^ + 7i)l log/ip. One can check easily that this is less than Ch~^^ 
unless u = 2 and /i < h~'^ ^^'^'^~^^\\ogh\^ in which special case there will be an extra 
I \ogh\^ factor (rather superficial); then we consider three-term approximation and one can 
see easily that the remainder estimate will be Ch~^'j. 

So we are left only with expression ()3.9|) for 

(3.10) u^u=-thY^ G'6{t)6{x2 - y2)S{xi - yi) 

?=± 

where here and below 'W means "replaced by" and two extra terms given by ()3.9j) again 
for 

(3.11) u^u^ = -ih G'A^G'5{t)5{x2 - y2)5{xi - yi), 

(3.12) u^U2 = -ihJ2 G'AiG'AiG'S{t)S(x2 - y2)Sixi - yi) 

?=± 

Actually factor | \ogh\ is superficial and the first term is of magnitude h~'^ji; further, all logarithmic 
factors below are superficial 
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where are corresponding parametrices and Ai = A — A. 

Also, note that replacing in Ai by Ai = (x2 — 2/2) ■ -Bi, -Bi = {9x2^)\x2=y2 o^i^ gets 

a proper remainder estimate with the exception of the special case in which in addition to 
terms due to ()3.10|) - ()3.12|) one must consider also an extra term given by ()3.9p for 

(3.13) u^^u'^ = -ihJ2 G'A2G'6{t)6{x2 - y2)S{xi - yi) 
with A2 = l{x2 - y2f ■ (5x2^) 1x2=^2- 

(ii) I will tackle terms due to ()3.11|) - ()3.13|) later; now let us concentrate on the main term 
due to ()3.10p . After this term is derived with T = Ch\ log/i| one can increase T to T\ where 
Ti = e/i7'^; however arguments of subsection 2.2 applied to A (whose coefficients does not 
depend on X2) imply that now one can take even T = +00. 

With T = +00 one can rewrite this term (before integration by ^/^ and integration over 
X2) as 

(3.14) (27r/i)~^ye(xi,xi;x2,^2,0,/i)(l -■?/'i(xi)) dxid^2 

where e{xi,yi]X2,C,2,^,fi) is the Schwartz kernel of 1-dimensional Schrodinger operator 
a(x2,'C2;^) obtained from A by change of variables Xi h-* 7x1 (and thus fi ^ 1) and 
hD2 ^ ^2 while e^{xi] X2,^2, K) would be Weyl approximation for its restriction to diagonal. 

However I would like to refer to a canonical operator ao rather than to more general 
operator a and thus few extra steps are needed. I will delay them slightly. 

(ii) Consider now the second term corresponding to 1)3.111) : I remind it does not exceed 
C/t7| log Zip with superficial logarithmiic factor factor. I claim that 

(3.15) The total contribution of the second term is modulo C fx^^^^h^^ + | logh\^. 

Really, note first that one can assume that fi < h~'^\\ogh\~^ , otherwise the total con- 
tribution C/i~^7i| log/ip is less than announced. Let us apply rewrite the term in question 
as 

(3.16) u^ui = -thJ2 G'BiG'[A, X2 - y2]G'S{t)6{x2 - y2)S{xi - yi) 

?=± 

-ihY, G'[x2 - 2/2, B^]G'5{t)5{x2 - y2)5{x^ - yi). 

Let us apply to ()3.16|) the method of successive approximations with unperturbed operator 
frozen as Xi = yi, then each next term will acquire factor C/i7'^-^/i| log /i| x (7/7i)''-i| log/i| 
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and therefore the second term does not exceed Ch~^'y\ log/i|~'^^ as long as 7 < 7I \ogh\^^. 
On the other hand the leading term of is just a Weyl expression without any auxiliary 
operators; this expression is odd with respect to (.^2 ~ Atx^'/z/) produces after integration 
with respect to ^2- 

On the other hand, as 7 > C7I log/i|^^ one can apply successive approximation method 
with unperturbed operator A is obtained from A by fixiu/u ^ C,2 + fJ'Zi~^{xi — zi), z\ = 
{fi~^h'^2y^'^ and the same is true for commutator [A, X2 — 1/2] and [Bi,X2 — 2/2]- Then 
since both xi and yi are confined to interval of the length (3 = (/i7'^~^)~^), each next 
term acquires factor (3'~f~^ \ log/i| x fi~^h\ logh\^~'' and therefore the second term does not 
exceed Cfi~^h~^ \ log/ip7^~'^ and its total contribution does not exceed what is announced 
as 7 > 7I log h\^^. 

On the other hand the leading term corresponds to harmonic oscillator and obviously 
sums to as Xi = i/i runs through TZ. 

So I will keep T = Ch\ logh\ for a while. In the analysis below I am going to consider 
fl3.9|) with r replaced by T' which contains integration with respect to Xi but not X2; then 
one can assume with no loss of the generality that 

(3.17) V, a, (j) do not depend on X2- 

It is convenient to have 0=1; one can always reach it now changing xi xi + 0{x1) 
and paying the price: now 

(3.18) A = ^(^a,ix,f{hD,f + a{x,f{hD2 - fix^u)' - V{x,)) 

with cr = cTi = 1 as Xi = 0. 

Let us apply the method of successive approximations again this time using as unper- 
turbed operator Aq obtained from A by freezing Xi = in V, a and ai: 

(3.19) Ao = ^ {{hD,f + {hD2 - ^ix'i/uf - y (0)) . 

In these successive approximations the first term is still 0(/i~^7i| log/i|) and an estimate for 
each subsequent term is gaining an extra factor CTo^i/h x 71 1 log/i|; therefore (n + l)-th 
term does not exceed Ch~'^'j^\ log/i|" which is small for large n since 71 < due 
to our assumption /i > h~'^ /(2»^-i)_ Note that we can write 

K 

with Ak = Ak{xi,X2, hDi, hD2; h, fi) as k > 1, 

k=0 
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Aq is already defined by ()3.19j) . Plugging it into successive approximation formula we get 
many terms but the first term results in 

(3.20) {27rhy^ J eo(xi, a;i; Xs, 6, 0, /i) (l - ^/'i(xi)) rfxid^2 

(after finally I take T = +oo instead of T = Ch \ log /i|) and the second term will be one for 
the approximation A = Aq + xiAi. 

But then the contribution to ()3.14|1 of this refined second term would be because it is 
equal to (before integration over ip dx2) 

(3.21) j Tri(6(6)V'i)rf6 

where ^('^2) is an operator in EI = L^(M), Tri is the trace in this auxiliary space and 

rb{^2)T = -biQ with = (-1)^6, (^tO(^i) = H-^i)- 

(iii) Now let us consider terms in this successive approximation with n > 3; they will be still 
too large just to be skipped. Let us apply to them the successive approximations method 
taking as unperturbed operator Aq with Xi/u replaced by y'(/u + yi~^{x — xi): 

(3.22) io = ^ [{hD,)' + {hD2 - fi{y'(/u + y^'ix - x,)))' - ¥(0)) . 

One can see easily that each successive term of this approximation is gaining factor fi'-f'^~'^T^ / h 
fi'~f'^~'^h'^\ \ogh\^. Further, contributions of the second terms to ()3.14j) will be equal again 
after easy calculations. Finally the remainder will be then less than 

Ch-^f^ X [fih^^-^-yi log/i|^ X C71I log/iK^ < 1 

where the first factor is the estimate of the 3-rd term in the previous approximations and 
the second factor is what we gained now. So we can leave only the first terms in these last 
approximat ions . 

Now when we sum all these terms with n > 3, we get (^S^^^ — S^"^^ (multiplied by 
"0(1 — '01 ) and integrated in the end of the day) where S^"^ is what we would get for 
the first term if we ran last round of successive approximations for it; and this is exactly 

S^^{xi,X2,T,^, h). 

However since we did not do it we have the same expression but with a correction term 



(3.23) (27r/i)-i j eo(xi/7; 6, r, h) - ^^o'''^ 



(xi,a;2,r,/i, /i). 
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(iv) Consider now the second term given by with u replaced by ui defined by (jH.llI) . 
This term does not exceed Ch~^'j\ log hl"^ (with a superficial logarithmic factor). Again 
applying the successive approximations of (ii) one can see that replacing G'' by (5q brings 
the extra factor 7|log/i|^ and therefore all the terms but the first one do not exceed 
Ch~^'yl \ logh\^ which is less than the right-hand expression of of ()3.5|) (unless in frames of 
the special case). Similarly one can replace Ai by 

Ai = -/i(x2 - y2)P{y2)xi'^^ {hD2 - fix'^/u) - ^(x2 - y2)w{y2). 
Or even let us consider the original second term; I remind it does not exceed Ch'^^'-)\ \ogh\. 

Let 

(v) Now there is a special case when one needs to consider extra terms due to (j3.12|) , ()3.13p 
and the second term in approximation of (iv); all these terms do not exceed C/i~^7^| \ogh\^ 
and one needs only to ret rid off this logarithmic factor. 

Replacing (5^, Gq by Gq linked to Aq defined by (|3.22p brings errors less than the 
right-hand expression of ()3.5|) . On the other hand calculations after this substitution show 
that these terms do not exceed Ch~^'jl (so there is no logarithmic factor) which even in the 
special case does not exceed Ch~^'^. Easy but tedious calculations I leave to the reader. □ 

Proposition 3.4. (i) Let condition (\1.21\ he fulfilled and let tpi be an even function, ip = 
V'(x2). Then as Ch'^ < < /i-'^V{2>^-i) 

(3.24) \{27ih)-' J eo{xi,xr,X2,^2,0,h){l-^Piix^))dx^d^2 

-J £^'^{x,,X2,C2,0,h){l - M^,)) dx,\ < C^h-\ 

(ii) Therefore estimate 5|) holds for Ch~^ < /i < /^-'^^/(^''-i) well. 

Proof, (i) instantly follows from proposition 13 . 21 for A = Aq and (ii) follows from proposition 
13.21 and statement (i). □ 

3.2 Inner zone. I 

Now I need to consider an inner zone ^inn = {l^^il < 7}. In this subsection I consider 

(3.25) h-' [ (F,^,-i,XTr(VV'i^))rfr 

J —00 
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with T = Ch\\ogh\ even if this formula does not represent r{ipipie)\r=o with a desired 
precision, leaving the necessary correction for the rest of the section. In the end of this 
subsection I will list cases (due to section 2) when provides the desired precision. 

Running an absolutely standard method of successive approximations (i.e. taking as 
unperturbed operator with x = y) one can recover Weyl asymptotics with C^h~'^ x = 
6*7"^ error with h= h/^ while the contribution of the whole inner zone to the asymptotics 
is of magnitude /i~^7 This 6*7"^ error does not exceed Ch~^^ as long as 7 > eh2 which 
means exactly that fi < Ch~2 . 

This is not extremely strong condition, much weaker than condition 7 > {h\logh\)3 
arising as h~^{h\\ogh\)2 < h'^'y to provide remainder estimate Ch~^'y for the spectral 
asymptotics in the general case. Moreover, it is even weaker a bit than condition 7 > 
{h\logh\)2 arising as h~^{h\logh\) < h^^^i to provide remainder estimate Ch~^'y for the 
spectral asymptotics under condition (j2.1(J3p . 

And in the inner zone I can write a standard Weyl approximation or magnetic Weyl 
approximation to my sole discretion^^^ 

Note that both Weyl expressions are actually more precise in the inner zone than I 
have written before, giving an error not exceeding C/i~^7 x h'^. Really, when we write the 
complete decomposition in powers of /i and h fixing 7 x fl~^^'^ but considering /i < /i as 
variable, we realize that terms with /x^^/i"^"'"^'^ all are in magnetic Weyl approximation, 
terms without fi are produced by Weyl expression for non-singular operator and thus are 
estimated by C/i~^+^'^7 rather than by Ch~'^h^^'^ and unaccounted terms contain factor 
^2fe^2« -^j^j-^ I > k>l and thus an error term does not exceed Ch^'^h'^'y = Ch?'^~^. 

Note that Ch'^'j"^ < Ch~^'j as 7 > h^^^ or equivalently /i < Ch"'^^/'^. Further, we get 



which can be rewritten as ()3.26p with 0(1) error. Therefore we arrive to 

Proposition 3.5. ^'^^ With an error Ch'^'j~^ one can rewrite (jj. 2^ with T = Ch\ \ogh\ as 




(3.26) 




Therefore due to subsection 2.3-2.5 we arrive to 



^^^It is easy to get rid off superficial logaritlimic factors in this case. 

i6)\Yhich would not be the case in the outer zone where we must use magnetic Weyl at least as ^ > h 
I'^^Cf. proposition El 
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Corollary 3.6. (i) In the general case the contribution of the inner zone to the spectral re- 
mainder estimate with the principal part ^3.2(^ does not exceed Ch~^'y + Ch~^{h\ \ogh\/'y)2 ; 
in particular as jj < {h\ log/i|)~§ the remainder estimate is Ch'^'j; 

(a) Under conditions . 1^2. lU!^ the contribution of the inner zone to the remainder 

estimate with the principal part !^S.2f^ does not exceed Ch~^'j + C\ log/i|7~^; in particular 
as II < {h\ Iog/i|)~2 this remainder estimate is Ch~^^; 

(Hi) Under conditions ^1.21\ . [ 2.91 ] with C, > Ch\\ogh\/'^ + C'^ the contribution of the 
inner zone to the remainder estimate with the principal part (j,3'.i:^6|) does not exceed Ch~^'y + 
C/i^7~^; in particular as fi < h~~i'^ this remainder estimate is Ch~^^. 

Then combining with the result of the outer zone analysis we arrive to 

Theorem 3.7. Let condition ^2.^ be fulfilled. Then 

(i) In the general case TZ defined by {\0.12^ does not exceed Ch~^^ + Ch~^{h\ log/i|/7)5; in 
particular TZ does not exceed Ch~^'y as fi < C{h\ loghl)^'^^^ ; 

(a) Under conditions 1^1.21^ and ^2.10S^ TZ does not exceed C/i~^7 + C7| \ogh\; in particular 
TZ does not exceed Ch~^^ as fx < C{h\ log 

(Hi) Under conditions and p.^7| ) with ( > Ch\\ogh\/^ + C7 TZ does not exceed 

Ch-^j + Cinh + l)V('^-i); m particular TZ does not exceed Ch~^^ as jj < Ch-''^/^'^''-^\ 

On the other hand repeating (almost all) arguments of the proof of proposition 13.31 one 
can prove easily 

Proposition 3.8. In frames of proposition 

(3.27) \h-^ I [Ft^h~^,XT{t)T{ijiJiu))dT- 

J —00 

'^^^(x, 0) - Sf^{x, 0))7A(x2)^i(xi) dx- 

(27r/i)"^y" ^x(-^^)) j Go{xi, Xi;x2, ^2, r, h)ip{x2)i'{xi) dxd^2 dr] < 

Ch~^^ + C\loghf 
Combining with the results of the previous subsection we arrive to 

Theorem 3.9. Let conditions (f7~^ and ( fO^ with ( > Ch\ log/i|/7 + C7 be fulfilled, 
= V'(X2). Then for /^-^'/(Si^-i) < ^ < h-"] \ogh\-^ 

(3.28) I I [e{x,x,Q)-£^^{x,Q)^i){x)dx- I {x2,Q)i){x2) dx2\ < Ch-^^ + C\\oghf 
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with 



(3.29) £^^ = i2nh)-' J eo(xi,Xi;X2,6,r,n)dxirf6- J Sf^{x,r)dx^. 
3.3 Inner zone. II 

In section 2 we derived asymptotics with the principal part which is the sum of 

(3.30) 5^/1-1 f U,^h-wXTjt)r{4j,tPQ^u)) dr 

where Qm are elements of partition (in hD2 and may be X2) and = Ch\ \ogh\/ g^. We 
can rewrite ()3.3U|) as the sum of ()3.25p with T = T = e^j and 

(3.31) Y.^'' i (^t^hMxTjt) - Xf(t))r(^i^Q™n)) dr. 

Consider now ()3.25p with T = e^y. Also due to propagation results of section 2 one can 
replace there T = e'j hj T = Ch\ log/i|; then all the results of the previous subsection are 
applicable. It includes ()3.27|1 as well where one can replace back T = Ch\ logh\ by T = €7 
in the first and/or the last term of the left-hand expression. 

Moreover again due to propagation results of section 2, all terms in fl3.31|) with > Cpi 
are negligible and therefore one can rewrite (modulo negligible) expression ()3.30|) as the sum 
of ()3.25|) with T = Ch\ \ogh\ and ()3.31|) with T = e7 with summation with respect to m 
such that Qm ^ Cpi- 

Proposition 3.10. (i) With an error not exceeding Ch^^'j + C7f ^ one can rewrite (jj. 3(J^ 
as the sum of i]3.26} and Jij) with T = e^y and Qm partition of unity in the periodic zone 
only. 

(a) In particular, as fi error does not exceed Ch 7. 

As /i > C/i-'^'/(2^-i) I need a bit more complicated analysis: 

Proposition 3.11. (i) With an error not exceeding Ch~^'y + C\logh\^ one can rewrite 
30^ as the sum of the following three expressions: 



(3.32) j (^S^'^ix^O) - £^'^ix,0)'^ij{x2)Mxi) dx+ 

J eo(xi, xi, X2, 6, 0, ^)^(x2)^(xi) dxdi2 
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where the first term is actually less than the remainder estimate, 

(3.33) {27rhy^ J (^^(~^^)) ^ '^(^)) J eo{xi, Xi;x2, ^2, r, h)ip{x2)'ip{xi) dxd^2dT 

with T = T = e7, Dirac 6-function and where now {Qm} is a partition of unity in 

the periodic zone only; 

(a) Further, one can rewrite 3^ as 

(3.34) J2^2nh)-' I {^H-lr)) - 6{r)) x 

eo(xi, xi; X2, 6, t, K}ilj{x2)ilj{xi)Qm{x2, 6) dxd^2 dr 



where now {Qm} is a partition of unity in the periodic zone only. 

Proof. Proposition follows from proposition 13.81 and propagation results of section 2 with 
the only exception that in ()3.34j) summation is taken over all m, and one needs to prove 
that the corresponding terms are negligible in non-periodic zone where condition 

(3.35) {\^2-Vhl\ + \d,,V\)>^p<Cp, 
is violated. 

Note that the expression in question is the sum of 

(3.36) {27rTy^x{-—T) / eQ{xi, Xi, X2, r, h)ip{x2)ipixi)Qmix2, C.2) dxd^2dT = 

j j (^Ft^h-wXTit){T'tpi{xi)Qmiy2,hD2)uo)^ijiy2)dy2 

with T running from 67 to +00 and Uq defined for operator Aq. 

If IC2 — ^^^^1 > p on partition element Qm then (j3.36|) does not exceed Ch^{T + 1)"** 
as T > Co/i| \ogh\/ just due to propagation results of section 2 applied to uq and Aq and 
since for CqHI \ogh\/ < ^7 as p > pi, these elements are covered. 

On the other hand if l^^^jV^I x p on the partition element in question, one can notice 
that Aq depends on y2 via —\W{ii2) only and introducing new variable y'^ = — |W^(|/2) we 
get mollification with respect to the spectral parameter and thus (|3.36|) does not exceed 
C/i*(T + l)"'^ as well. I leave details to the reader. □ 
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3.4 Periodic and near-periodic orbits. I. Pilot-model 

Thus I need to analyze periodic zone Zp^r more accurately defining it by ()3.35|) rather than 
by 

(3.37) \^2-Vhl\^p<Cp, 

in the case when magnetic field is strong enough to prevent corollary 13 . 61 from presenting a 
sharp remainder estimate. 

Namely I need to consider term 

(3.38) h-' I (Ft^h-^r{xT{t) - Xf (t))r(^i^gM)) dr 

with partition element Q = Q{x2, hD2), T = Ti where Tq and Ti are defined in subsections 
2.8-2.9 and T = CqHI log/i|. Here due to propagation results of these subsections I can take 
any T G [Tq, Ti] to my discretion, in particular T = Tq. 

Let consider a pilot-model first; more precisely let us assume that operator under con- 
sideration coincides with a pilot-model in -6(0, 1). Then Tq = Cp~^h\ logh\ and Ti = ep 
with 1^2 — P on suppQ. 

Then fl3.38|l will be the same (modulo negligible) if one replaces m by m constructed for 
operator coinciding with pilot-model in M^. But then one can replace T = ep~^ by T = oo. 
Also one can remove restriction from below on p : p > 7. 

Also in this case one can replace ipi by 1 since contribution to (j3.38|) of the outer zone 
is negligible due to propagation results there and contribution of the forbidden zone is 
negligible as well. Further, one does not need to integrate over %l){x2)dx2 in (j3.38|) anymore. 

After all these modifications there is no need to consider a partition and then correction 
fl3.38|) is transformed to 

(3.39) h-^ I (Fj_^-v(l-Xf(i))rM)rfr 

J —00 

where I remind Vv = J v{xi, X2] Xi, X2) dxi, in our case does not depend on X2. 
Obviously 

(3.40) h-' (F,^,-i.rn) dr = (2nh)-^ j no(6, K) d^2 

with no (.^2; ^) the number of negative eigenvalues of operator a = ao(.^2; ^) with W = 1. 
On the other hand, 

(3.41) h-' r (t)r'M) dT= f S^'^ix, 0) dxi mod 0(71"') 
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and therefore 



(3.42) 



MW _ 



(27r/i)-i / no{^2,h)d^2- / ^o^(x,0)rfxi 



corr 



The error in ()3.42|) does not exceed Ch'^'j even if 7^^^ > h~^'y because in this case I aheady 
made an error replacing left-hand expression of ()3.41|) by its right-hand-expression in the 
analysis of subsection 3.1 and now I just compensated it by adding the skipped term. One 
can consider ()3.42|1 as a definition of the correction term for the model operator. Thus we 
arrive to 

Proposition 3.12. For operator coinciding with pilot model in (—1, 1) x M 



in (— |, |) where is defined by l \3.4^ - 

Expression (j3.42|) . multiplied by ip{x2)dx2 and integrated will be used in two next subsec- 
tions in the general case as well. However to get more explicit even if less precise expression 
one can replace ()3.39|) by 



(3.43) 





(3.44) 




with (p supported in [— 2Co, 2Co] and equal 1 in [—Co, C — 0], making 0(^1 ^) 
Using 



error. 



(3.46) h-' / (Fi_,-i.xf r'0(/il^2)w) dr = {2TihY^ / 71^(^2, ^)0(6) 



(3.45) / {Ft^h-^,T'cP{hD2)u) dr = {27rhr' / no(6, /i)0(6) d^2, 




—00 



—00 



with the second equality modulo 0(7 ^/i^) with 




the Weyl approximation of no(^2, ^), we arrive to 
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Proposition 3.13. For operator coinciding with pilot model in (—1, 1) x M 
(3.48) = (27r/i)"^ I (no(e2, h) - n^i^,, h))<P{^,) rf^ mod 0{-f^') 



in (-ii). 



To calculate ()3.48|) more explicitly I need to calculate eigenvalues A„(.^2, ^) of operator 
^o{C,2, K) and I am interested in those eigenvalues which are close to as .^2 is close to /c^. 
These eigenvalues are defined modulo 0{h^) from Bohr-Sommerfeld condition 

(3.49) i'^^^^, r) + ^-iM (/:(6, r)) G Z as r = A„ 

where 



(3.50) S{i2A)=i iidx 



and lm{'C.{^2, t)) is Maslov' index of the closed trajectory C{^2, t) (on 2-dimensional phase 
plane) on the energy level r; iM{'^{C2, t)) = 2 as r ^ 0. Since drS{^2, t)| = ^"(^2, 't) due to 
Hamiltonian mechanics, the spacing between two consecutive eigenvalues is 2'Kh/T{^2,'T) + 
0{h^). Further, all these eigenvalues A„(^2, K) are uniformly analytic functions of ^2, |^2| ^ 
1, <^ 1 as |A„| ^ 1 and also d'^^Xn > eo. This yields immediately 

Proposition 3.14. As ^ < e{h\\ogh\)~'' 

(i) iCr^l < Ch-^h-2; m particular, \Ef^^ \ < Ch~^^ as fi < Ch^'^/'^; 

(a) One can calculate with an error not exceeding Ch^^h (which does not exceed 

Ch"^'j as jj < Ch""^'^) replacing An(^2,^) by its approximation modulo 0(r]^ + fi?) where 
here and below f] = C,2 — kn! 

(Hi) One can calculate with an error not exceeding Ch~^h^^'^ (which does not exceed 
C{h~^'y + 7f^) as u >?)) replacing \n{^2, ^) by its approximation modulo Oi^rf + h?); 
(iv) As V = 2 one can calculate with an error not exceeding Ch~^h^ (which does not 
exceed C'ji'^) ) replacing Xn{^2, ^) by its approximation modulo 0{rf + h^) . 

To exploit (ii) one can replace A„(^2; K) by its approximate value 
(3.51) A„(e2,^) = A(6) + ^^^^^^. 
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Here A (.^2) is a "classical eigenvalue" corresponding to n = and Maslov index (instead 
of 2) and one can find it from ()H.49j) as r = 0; then A(^2) = — 'S'(.^2, ''")/^r>S'|^^g. 

I remind that {drS)\^^^ = T(^2), %5 = with 1(^2) defined by for even 1/ 

and similar formula for odd u; further, I{k*) = 0, d^^I{k*) = k > 0. Then 



(3.52) 



£. 



MW 



:(27r/i)-i 
{2nh)-^ 



2-Kh 
1 

2Txh 



.2Txh 



di2 



{So + i:^r^') 



2Txh 
mod 0{h~^h) 



drj 



with function G defined by 
(3.53) G{t) 



1 2 1 



with the converging integral in the right-hand expression. 
One can prove easily that 



(3.54) 



G^O, G(t+l)=G(t), ! G{t)dt = 0, GeCi 

Jo 



(3.55) 



To exploit proposition 13. 14( 111) . fiv) one can use more precise version of (j3.51|) 

{2n+l)7rh 



A„(e2,/i) = A(6) 



^(6 



and after obvious calculations ()3.52p is adjusted to 

(3.52)* = i27rrlh~'h^K'^G(p- + K,h) 

(3.52)** 8^^ = {27ryh-^ h^K-'^ (^G{t) + ushGit) + K^hG^it) 



mod 0{h~^h^-)., 



27ih 

mod 0{h~^h^) 



with some constants Ki, . . . , ^4 and 



(3.56) 



Gi(t) 



G(t') dt' 



t')G{t') dt', 



Gi E C'^/^ satisfying (j3.54p . Thus we arrive to 
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Proposition 3.15. For operator coinciding in B{0, 1) with pilot-model with potential V = 1 
and fixed tp, ipi (tpi = 1 as \xi\ < e) in asymptotics ( |J.^7| ) 

(i) Modulo 0{h-^h^^^), 0{h~^h), 0{h~^h'^ '^) andO{h~^h^) one can respectively skip cor- 
rection term 8^^ or to define it by (TOl . |(3.52)*[ |(3.52)**( - 

(a) In particular, uncorrected asymptotics has sharp remainder estimate if and only if 
m <C^ I.e. fi < Ch-"/^. Further, one can skip define by without spoiling 

remainder estimate as Ch~'^^^ < < Ch"'^^'^. 

Remark 3.16. (i) One can then generalize proposition 13.151 immediately to a pilot-model 
with potential V = W = const instead of = 1 by replacing h i— > hW~^^'^, jj, ^ fiW'^^'^ 
and therefore h i— hW~^'^'^^^^'^'^ just modifying (j3.52p to 

(3.57) - (2.)-ift"S*K-4»yS?G(^^) 
and similarly modifying |(3.52)*[ |(3.52)**| 

(ii) One should expect the similar correction for the general operators as well but one can 
see easily that under condition l^j^^iyi > e integrated correction term 

(3.58) J 8^^{x)ijdx2 

is negligible, and under condition ()2.103p it should not exceed Ch~^h (and be of this 
magnitude if there are critical points of W) with integrated versions of dSSai, 1(3.52)*! 
valid modulo 0{h~^h2), 0{h~^h?) respectively. 



3.5 Periodic orbits. II. General settings 

Now let us consider inner zone in frames of subsections 2.8-2.9 when we managed to prove 
sharp remainder estimate C'jh'^ as /i < h^~'^ producing the final result as the sum of 
expressions ()3.3()|1 taken over partitions ip = i^m{x2) and Q = Qm{hD2) of unity; I remind 
that ipm is -^-admissible and Qm is p-admissible with £, p satisfying some conditions the 
most important of which are ()2.70|) - ()2.72|) : here i denotes the function associated with the 
finest of subpartitions and we assume that on this element 

(3.59) \WV\ < CoC, 16 - klV{'d,X2)\ < C^p 

with one of these inequalities being reversible (with eo instead of Co); then one can take in 
()3.30|) any T ranging from 

(3.60) To = Ch\ log/i| min(4, ^) = CC^h\ log/i| x Ti 
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to Ti. 

There are also exceptional elements with T = T = Ch \ \ogh\ but they will be treated 
easily in the same manner. 

Therefore I need to consider expression ()3.30|) with T = T = Ch \ logh\ and also correc- 
tion terms ()3.38|) with arbitrarily T e [To,Ti] depending on the partition element; index m 
indicating partition element I am skipping. 

One can replace then [xrit) — Xfit)) by the sum of Xrit) with T running from T to Tq 
and with x supported in(— 1,— |)U(|,1) thus replacing ()3.38p with the sum of terms 

(3.61) h"' I (Ft^h-^rXT{t)T{iPiijQu)\ dr = iT'^ Ft^h-^rXT{t)T{il)i^Qu) 

J — oo 

with x{t) = t~^x{t) and T G [f , To]. 

Further, one needs to consider only T ranging from T' = e7 to Tq because terms ()3.6H) 
with T ranging from T to T" are negligible; therefore partition elements with Tq < T' 
will be eliminated. 

Remark 3.17. (i) Obviously (|3.61|) does not exceed C/i~^7p£ and thus (|3.38|) does not 
exceed Ch~'^'ypi\ log/i|; 

(ii) Further, one can decompose ()3.61|) in the same way as in the proof of proposition 12.301 
therefore (|3.6H) does not exceed Ch~'^^pl x h\ \ogh\/^ = Ch~^pi x | log/i| and thus (|3.38p 
does not exceed Ch~^pi\ log/ip. 

Even if this is the remainder estimate obtained by taking T = e7 in the Tauberian 
arguments (with a superficial logarithmic factor), it provides a solid ground for the further 
estimates. 

Let us consider regular elements and apply successive approximation method freezing 
in unperturbed operator X2 = y2 as usual. However to estimate approximation terms I will 
not commute (x2 — ^2) with G^, but simply remember that \x2 — y2\ ^ CpTo + C7 on 
the partition element in question in propagation as T < Tq; this is due to proposition 12.241 
extended according to proposition 12.331 

In this subsection I consider the case when condition ^2.1U'J^ is fulfilled leaving the most 
general case for the next one. 

Then p = ( = i and Tq = Cp~'^h\ \ogh\ and hence \x2 — 2/2! ^ CpT^ + 7. Therefore 
perturbation Ri = (x2 — y2)dx2^ + 0{\x2 — y2p) does not exceed 

CC(pTo + 7) + C(pTo + 7)2 X r ""^ {h\ \ogh\ + p7) + P"'/^'l log/ij^ + f ). 
Or cquivalently with (C^ + p^)^ > Ch\ \ogh\. 
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Therefore while the first term of what is obtained when one plugs approximations into 
()3.61|1 does not exceed Ch~^p£\ logh\ = Ch^^p'^ \ \ogh\, each next term acquires a factor 

(3.62) ryiT/h^ {h\\ogh\ + p^ + p-'^h^\\ogh\'^ + f) x p~'^\\ogh\] 

to keep this factor less than 1 one needs to consider 



(3.63) p> pI = Chi\\ogh\ 

which is greater than C7 because under condition ()2.1()3j) and 7 > Po remainder estimate 
0{h~^'y) is already proven without any correction terms. Under condition (j3.63|) expression 
fl3.62|) is of magnitude of p~'^h\ log/ip. 

Now the second term of approximations does not exceed 

Ch~^p^\\ogh\ X p-2/i|log/i|2 X \ \ogh\^. 

Finally, summation with respect to all partition results in C \ \ogh\^ which does not exceed 
Ch'^'y under condition ()2.11ip . 

On the other hand, contribution of core {p < p^} to the remainder estimate due to 
Tauberian arguments with T = e'j does not exceed Ch'^pl"^ < C\ \ogh\^ . 

So, one needs to consider only the first term of approximation. Therefore, under condi- 
tions ()2.103p and p < h~'^\ \ogh\~^ remainder estimate Ch~^'y still holds with u replaced 
in the principal part ()3.30|) by the first term of this successive approximation procedure. 

Thus we can take Tq = 00 thus arriving to 

(3.64) = y"((27r/i)-ie(xi,xi,0;x2,6,^,/^)- 

/ Ft-,h-^rXf{t)'^x{U) dTj'i/j{x2)'ipi{x)(p{^2)dxidx2d^2 
J —00 

where e = e(xi, yi, r; X2, /i, /i) is a Schwartz kernel of the spectral projector, associated 
with one-dimensional Schrodinger operator a = a(x2,^2,^, p) which is obtained from A 
by replacement hD2 ^ ^2 and subsequent change of variables xi 1— xi/7 but not setting 
a = (f) = l, V = W and U = U{xi,yi,t;x2,C2,Kf^) is a corresponding propagator, i.e. 
Schwartz kernel of e^^ Here and below Q = ip{x2)^p{C,2) denotes partition element. 

However after summation over xi-partition remembering that as > C on supp'i/'i 
we get negligible terms, we can replace ipi by 1 thus arriving to 

(3.65) £^^^Q = j(^i2nhr'ni0;x2,^2,h,p)- 

h-' f Ft^t,-.,Xf{t)T'{U)dT)^{x2)^{i2)dX2di2 
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with n eigenvalue counting function associated with the same operator a; I remind that T' 
includes integration with respect to xi. 

Without the last term formula ()3.64|) delivers the correct principal part of asymptotics 
which after summation over all partitions (including non-periodic and outer zones) is 

(3.66) (27r/i)~^y e{xi,Xi,0;x2,^2,fi,IJ')'4'{x2)i'iixi)dxidx2d^2 

which as we show in the next section is a correct answer without condition (j2.1(J3j) as well. 

However, I would like to replace the reference to operator a by the reference to operator 
ao. To do this I need the following 

Proposition 3.18. Let condition 1^2. lU!^ be fulfilled. Let p > Cp^. Then 
(i) Eigenvalues A„ o/a satisfy inequality 

(3.67) \^x2,i2^n\ ^ P 

as \\n\ < ep^ ; 
(li) As pyCj 

(3.68) |A„(x2,6) - A°(x2,6)| < is even, 

(3.69) |A„(X2,6) - A°(X2,6) - /5(x2,e2)6l < Cf , \(3\ < vis odd, 

where A° are corresponding eigenvalues of ao . 

Proof. Note that a = ag + ta' with a' bounded by C7(a + Cq) and t = 1 therefore 
due to Rellich decomposition one needs to prove that |9jA^| < C^'^ as t = 0. As is 
even eigenfunctions T„ of ao are either even or odd and therefore (a[jjjT„, T„) = while 
|«_T„,T„)|<C7^ 

Consider odd v. Then /c^ = just due to symmetry and the above arguments are appli- 
cable for /c = 0. On the other hand A„ = A„^sc, A° = A° g^. mod 0(/i^) with semiclassical 
expressions A„,sc, A,° g^. and we need to prove ()3.69|) for these expressions; due to analyticity 
etc these equalities are true as 7 > Ch\\ogh\ i.e. in our assumptions. Since A^^^, is ana- 
lytic with respect to h,t,^2, coincides with A° g^, modulo 0(7) and also coincides with A° 
modulo 0(7^) as ^2 = 0, we get (jSSHl)- ' '□ 

I want to remind that also 

(3.70) C'^Q = /( f Ft^hMxrit) - Xf W)r'(t/) rfr)^(e2)V^(x2) ^^^2^6 
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and one can rewrite it as the sum of 

(3.71) j (y"° Ft^,,-i,m^'iU)dT)v{^2)H^2)dx2d^2. 

Proposition 3.19. In frames of proposition W~TR as h""/"^ < /i < h~'^ \ \ogh\~^ 

(i) With an error not exceeding Cfrfh~^ one can replace in the right-hand expressions of 
< \3. 7(J^ tJ by Uo associated with operator ao(x2,C2, therefore in (),'/. 6',')|1 one can replace n 
by no and simultaneously U by Uq; 

(a) < 6*7"^ and this estimate cannot be improved unless V{0,X2) has no critical 

points. In particular, as ji < Ch~'^^'^ and only then one can skip without deteriorating 
remainder estimate C'^h~^. 

Proof. I remind that under condition (j2.1(J3|) the contribution of the partition element to the 
correction term does not exceed C/i~^p^| \ogh\^ . Then running successive approximation 
method with unperturbed operator ao and using proposition 13.181 one can see easily that 
each next term gains factor Cp7 x Tq/H = Cp"^7| logh\ and thus the error does not exceed 

Ch'^p^l loghf X P-I7I \ogh\ = Ch-^p^l loghf 

and summation over partition results in C/i~^7pi| logh^^'^'^ ^ h~^^. 

Then (ii) follows from (i) and the calculations similar to those of the previous subsection. 
I leave easy details to the reader. □ 

Then we arrive to 

Proposition 3.20. Under conditions 1^2. lU^ as h'"^'^ < p < h'"] \ogh\~^ 
(i) Asymptotics 

(3.72) [ e{x,x,0)i!{x)dx= f £^]^ {x)i){x) dx + I {x2)i){x2) dx2 modO{h^^^) 



holds with defined by 

(ii) Moreover, with 0{h~^K^^'^ + 7f^)-error one can replace S^"^ by expression [ 3.51 ] and 



with 0(7^^ -error one can replace by expression \{^.b2)*\ modified in the same way. 
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3.6 Periodic orbits. III. General settings (continuation) 

Let us consider more general settings. Then, as I have shown in the previous subsection, for 
any regular partition element described in subsection 2.9 its contribution to the correction 
term does not exceed Ch~^ip\ \ogh\^ (where actually K should be equal to 1). 

Let us consider a successive approximation method with unperturbed operator obtained 
by fixing X2 = y2 in the coefficients. Again as before \x2 — 1/2 1 does not exceed C{pTq + 7) 
and perturbation is (x2 — 1/2)^ with B = {dx^A\x^=z2 calculated at some point Z2 of the 
same partition element (containing both y2 and X2) and thus can be estimated by C(; thus 
the perturbation is estimated by C({pTq + 7) and the contribution of the second term to 
the correction term is estimated by 



(3.73) Ch~^pi\ log/i|^ X C(pTo + 7) X To/h < 

cic[p\p' + ar' + h~'jp{p' + ar') \ iogh\ 

+ /i-i7p)|log/i|^+2 < Ch~^^logh\ 



\K+2 ^ 



K+3 



as To = C/i(p2 + (^)-i|log/i|. 

This was a contribution of the regular {£, p)-element. This element was actually a subele- 
ment of some more rough element and then the contribution of this larger element would 
be a sum over finer subpartition. So in frames of subsection 2.9 as long as we use nondegen- 
eracy condition |(2. 105)"^ contribution of this larger element would be Ch~^^2 \ \ogh\^~^'^ as 
well. Continue this process until the very top we get Ch~2^2 \ \ogh\'^^ with ^ ^ thus 
resulting in C\ \ogh\^. So, under nondegeneracy condition |(2.105)^;„| we arrive to a proper 
error estimate 0(/i~^7) as /i < h~'^ \ \ogh\~^'^. 

On the other hand, without nondegeneracy condition |(2. 105)"^ we still have this condi- 
tion fulfilled until the very top when we join regular elements to -8(0, 1); then we arrive to 
the error estimate 0(/i~^7 + h~^^ with arbitrarily small exponent S > 0. 

On the other hand contribution of the irregular elements to the correction terms are 
less than C\ \ogh\^, Ch~^ with/without non-degeneracy condition [(2.105)^ respectively. 

Thus we arrive to 

Proposition 3.21. Let p < Ch^"] \ogh\~^. Then 

(i) Under nondegeneracy condition \{2.1Qb\^ asymptotics holds with the principal part de- 
fined by 6'6|) and the remainder estimate 0(/i~^7); 

(a) Without nondegeneracy condition |(2.105)"^ asymptotics holds with the principal part 
defined by fIJ. 6'6|) and the remainder estimate 0{h~^'^ + h~^^ . 
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Our next step is 

Proposition 3.22. Let p > Cj. Then eigenvalues A„ of a satisfy inequalities 
(3.74) P~'|V5,A„|+C-'|V.,A„| <C 

and (fJ^ or (fJ^ as |A„| < ep^ 

Proof. Proof of (j3.74|) is obvious and (j3.68|) . (13.691) are already proven. □ 
Proposition 3.23. Let p < Ch'^l \ogh\~^. Then 

(i) Under nondegeneracy condition |(2.105)"^ asymptotics ( | J. 7^ holds with the correction 
term defined by l\3.4^ and the remainder estimate 0{h~^'^^; 

(a) Without nondegeneracy condition \(2. 105)^ asymptotics 72 holds with the correction 



term defined by l\3.4^^ and the remainder estimate 0{h -'^7 + h ^) . 
(Hi) Moreover, with 0[h~^h + ^y^^)- error one can replace by expression (jj. and 

with 0(/i~^/?.^/^ + 7f"'^)-, 0{yi^) -error one can replace by expressions \{3.52)*\ |(3.52)**| 
respectively modified in the same way . 

Proof. I remind that contribution of the regular (£, p) partition element to the correction 
term is 0{h~^ptj. For even u application of the successive approximation method brings 
in every step factor 6*7^ x T^/h = C7^p~^| log/i| and the contribution of the second term 
is 0(/i~^7^p~^£| \ogh\y Then summation over partition with p > po = C'tI log/il^ results 
in 0(/i-^72po^l log^l^) = 0{h-^^). 

For odd i> application of the successive approximation method brings in every step 
factor C^p~^\logh\ and all above estimates remain true for the third term in successive 
approximations and even for the the second term with the single exception of its part which 
is of magnitude 0(^h~^'yi\ \ogh[j and which emerges from the linear with respect to ^2 term 
in formula (j3.69p . One can see easily that this exceptional term becomes after integration 
with respect to ^2- 

Now in the both cases one needs to look at contribution of zone {|'C2 — < Po} which is 
0{h^^PQL) where L is the total length of the corresponding X2 intervals and thus is 0{h~^'^) 
under condition |(2.105)nl| otherwise it has a superficial logarithmic factor we need to get 
rid of. 

When looking at such narrow zone due to above arguments one can skip all the 0{p^) 
terms in A„ — A° and thus we can assume with no loss of the generality that both A^ 

Cf proposition EH 

Surely under condition 12. 1041 better estimates (depending on m) holds. 
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and A° are quadratic forms with respect io t] = ^2 — and 7 with coinciding coefficient 
at rf and after obvious transformations we can assume^^^ that A° = xrf + xq, A„ = 
xrf + Xi with Xi — xq = 0(7^). Then we can take integration with respect to .^2 all over 
M. But then both correction terms would be of the form cih~^h^^'^G{c2 + c^xj/h) with 
j = 1,2 and equal coefficients Cfc, k = 1,2, 3; since G E C^/^ the difference does not exceed 

ch-^ny^ X (7Vn)i/2 = Ch~^^. □ 



Thus all main theorems are proven as /i < /i log /i | 



4 Superstrong Magnetic Field 

In this section I assume that fi is close to h^'^; more precisely I assume ffist that 
(4.1) /i"^"^ < /i < e/i"^ 

with very small exponent S > constant e and I leave the case eh~'^ < < Coh'" for 
subsection 4.5. Under assumption ()4.H) parameter h = h/'y is small but not very small: 

< h <^ 1. Moreover, as h> e\ logh\~^ variable xi is no more microlocal. 

I am going to consider operator in frames of pdo theory with operator valued symbols 
in the axillary space H = L^(M) with inner product (., .) and norm |.|. 



4.1 Outer zone 

Propagation estimates in the outer zone are basically done: as /i < eh~'^\ log/ij"*^ I already 
proved that the propagation trace Ft^h-i^XTif)^{'^Qu) is negligible as T G [TojTi] with 
To = Ch\ \ogh\ and Ti = e'p, where ip = 'ip{x2) and Q = Q{hD2) with p < \^2\ < 2p on the 
support of Q and C < p < C[fih'^) ^^^'^ ^\ Further, as p > C^ph") ^^^^ the standard 
ellipticity arguments imply that the propagation trace is negligible even as Xt^^) is replaced 
by XT^t)- 

Also I proved that as eh~^\ log/ij^*^ < /i < Ch~^ propagation trace Ft^h-iT-XTit)r (ipQu) 
is negligible as T G [To,Ti] with the same Tq,Ti and C < p < C\logh\'^ while ellipticity 
arguments work as p > C| logh\'^. 

But then one can prove estimate 

(4.2) \Ft^h-wXT{t)r{i'Qu)\ < Ch-^-ip^/-\ \ogh\ 

which is "almost perfect" : the problem here lies only with the logarithmic factor. 
^^^This assumption is valid only if we arc interested in the difference of correction terms. 
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To improve estimate ()4.2j] let us launch launch successive approximation method, fixing 
in the coefficients of unperturbed operator X2 = y2 as usual and also fixing xi = there 
but not in fiXi/u. Then one can estimate easily the final contribution of the first term by 
the same Ch~^'jp^^'^ \ log/i| as before and also estimate by C contribution of all other terms. 
Therefore in the estimate part only the first term should be counted. But then one need to 
consider 

(4.3) h^' J (F,^,-i.(r'f/(t;6)))Q(6)rfa;2rf6 
where 

(4.4) U = U{xuyi, t- X2. 6) = e'""""("^'«^''^^*T„(xi; xa, 6, ^)T„(2/i; xs, 6, K) 

n 

is a propagator for 1-dimensional operator which after change of variable xi t— Xi/7 is just 
ao = ao (0:2,6,^). In the below arguments let us skip integration over X2 and X2 in nota- 
tions . So, A„,(6,/l) are just eigenvalues of ao(.^2,^) and T„,(xi; .^2, ^) are orthonormalized 
eigenfunctions. Plugging ()4.4|) into ()4.3p I arrive to 



(4.5) 



J^h^'T J ^(^(r-A„(6,/^)))Q(6)rf6 



and in the case of even u I will consider separately sums with respect to eigenvalues corre- 
sponding to eigenfunctions odd and even with respect to Xi. 

Proposition 4.1. As h <ti 1 eigenvalues Xni.^2, K) of ^^{^2, have the following properties 
as 161 > Ci.- 

(^) Xn{^2,h) > A"^ eoh\^2\^'^-'y'^ - Co; 

(a) For even v and 6 < — e all eigenvalues are larger than Co^; 

(Hi) Spacing between two consecutive eigenvalues'^'^^ \n{^2,K) and \n+i{^2,K) belonging to 
interval (— eo, eo) is not less than A; 
(iv) Also for these eigenvalues 

(4.6) Co>6%A„(6,^)>eo, 



As V is even I consider two series of eigenvalues separately: those with even with respect to Xi 
eigenfunctions and those with odd eigenfunctions. 
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Proof. Proof of (i)-(iii) follows from Bohr-Sommerfeld theory with the semi classical param- 
eter h' = hp^'^"^^^'' when it is small enough. Note that zone {|xi| < e^yp^^'^} is classically 
forbidden with the smaller parameter h" = h/ p. 

The easiest way to prove ()4.6p is to note that propagation results of subsection 2.2 
basically are equivalent to these inequalities as fi < €h^''\logh\~'^] however this condition 
could be replaced just hj h < e' because eo{xi,yi, ., .) depends on h, C,2 rather than /i, /i 
separately. □ 

Let us consider ()4.5|) . Since x is fast-decaying at infinity, each term in ()4.5|) does not 
exceed 

h-'T j{l + T\T-X4C2)\h-yd^2-p 

due to (jHH); on the other hand the number of terms in ()4.5p does not exceed CA~^ x 
;i-ip{i-^)/^ Therefore 

Proposition 4.2. As h <^ 1 and p > C expression does not exceed CfT^p^^^ . 

So, I got rid of the logarithmic factor in ()4.2j) . Then dividing by Ti = ep \ conclude 
that the contribution of this partition element Q to the remainder estimate does not exceed 
Ch~^p^^~'^^^'^. After integration over dp/p I conclude that the contribution of the whole 
outer zone to the remainder estimate does not exceed Ch~^ = Ch~^'y. So I arrive to 

Proposition 4.3. Under condition ( |^.iD contribution of the outer zone ^out = {\^2\ > C*} 
to the remainder estimate does not exceed Ch~^'y. 

4.2 Inner zone. I. Estimates 

So far I have no results in this zone Z^^^ = {\C,2\ < C'o} as /x > eh~'^\ \ogh\~'^. Thus I cannot 
even refer to operator ao directly. However let us start from the properties of ag = ao(^2, ^) 
with potential V = 1. Let us first to move through ^2 = 1 or ^2 = ±1 for even/odd u 
respectively. 

Proposition 4.4. Let h < e and 1 — e < |^2| < C'o with small enough constant e > 0. Then 
operator ao(^2, ^) is microhyperbolic with respect to ^2 on energy level 0; more precisely 

(4.7) ^2{d^M^2,h)v,v) > e\v\^ -Ch-'-\aio{^2.h)vf \Jv G /^(ao) 

with large enough exponent r. 
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Proof. Let us consider ,^2 > 1 — e in the both ,^2 < — 1 + e and odd u is treated 

due to the symmetry. It is sufficient to prove ()4.7j) for v G (E(,^2, ^; e^'') — E(,^2, ^; — efi''))EI 
where E(,^25 ^; A, /i) is the spectral projector of ao(.^25 Then the standard WKB approach 
imphes that for any v described above with = 1 and e > one can approximate v by 
WKB solution as Xi > e: \v — w\j^^^y^^^ < df'"^ where w is a. WKB solution constructed 
for A = 0, C = C(ei), and \w\^^^^^^y~< 2. Then 



\w\ = X — , — 0(1 

with < X < c. 

Then, from explicit the WKB calculations it follows that 

{d^,ao{^2,h)w,w}{^,>e} = const ( / == - o(l) j > 3e|w;|^ 

and therefore 

(4.8) (%ao(6,%,t;){.,>.}>2e|^|J,^>,^-Cr-2. 

The last inequality holds for domain {xi < — e} as well because it is either classically 
forbidden (for odd u) or just due to symmetry (for even u). Furthermore ()4.8|1 obviously 
holds for domain {|a;i| < e}. After summation I arrive to ()4.8p with integration over M 3 Xi. 
Easy details linked to WKB calculations I leave to the reader. □ 

Proposition 4.5. h < ho{eo). 
(a) Further, let v he even. Then 

(i) As 1 + e < ^2 < Co spacing between eigenvalues'^^'' Kii^, ^) o-nd A„+i(^2, K) (since v is 



even I separate two series of eigenvalues; see footnote IffZjj at least eh where here and 
below e is a small positive constant; 

(a) As eo < ^2 ^ ^ — ^0 spacing between eigenvalues A„(^2,^) o,nd Xn+i{^2,h) (I do not 
separate two series of eigenvalues anymore) is at least eh; 

(Hi) As 1^2! < Co spacing between eigenvalues Xni^2, h) and Xn+i{^2, h) is at least eh?^^^"^^'' ; 

(b) On the other hand, let v be odd. Then 

(iv) 74s 1 + e < 1,^2! ^ C'o spacing between eigenvalues Xn{^2,h) and Xn+i{C,2,h) is at least 
eh; 

(v) As \^2\ ^ ^ — ^0 spacing between eigenvalues A„(^2, h) and A„+i(.^2, h) is at least eh. 
^•^^ In this proposition I consider only eigenvalues belonging to [— e, e]. 
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Proof. All cases but (iii) follow easily from the standard WKB method. Moreover, (iii) 
follows from the standard WKB method and rescaling. □ 

So spacing between eigenvalues is rather large - larger than eh? (as h"^ > h~^) and 
therefore their derivatives of order a could be estimated by Cah~^°'. Further, \d^2,fi/^n{C,2)\ 
is bounded. 

Then decomposing u{x, y, t) into T„(xi/7; hD2, h) "^^^ which are functions of xi but also 
h-pdo theory with respect to X2 

(4.9) u{x,y,t) = y^/f~'^Tn( — ]hD^^,h)uni{x2,y2,t)Ti( — ]hDy^,hy 

one can rewrite the basic equation {hDf — A)u = as a system 

(4.10) (hDt- An{x2,hD2,h)yura = ^Bkiu^i, 

k 

(4.11) A^{X2, hD2, h) = X^{hD2, h) - \W{X2) + \ 

with operators Bki estimated by Chh"^ where I increase K if needed. 

Let us consider e^^-sized intervals with respect to X2 and r. On each such interval at 
most one of operators hDt — An{x2, hD2, h) fails to be elliptic there (with the ellipticity 
constant e/l^). There is one exception: if u is even and ^2 > 1 — e "odd" and "even" 
eigenvalues^^-* could be pretty close to one another and then there are at most two such 
operators with numbers n and n + 1 {I cannot separate eigenvalues anymore; n corresponds 
to even eigenfunctions and (n + 1) to odd ones). 

Let n be such exceptional number. Then one can rewrite system ()4.10|1 either as a single 
equation 

(4.12) (hDt - A^{x2, hD2, h) - B'^Uni = 

where Ukn with k ^ n are expressed via Uni and B'^ = Bnn modulo operators not exceeding 

An ^ 

. A„+i^ 
operator B'^. 

Proposition 4.6. In frames of proposition ±d^2^n{^2, h) > ei as long as ±{C,2 — k*) > 
eo and |A„(^2)| < ei- 

^'^^ Where T„(a;i,^2, fi) are eigenfunctions of ao($2, ft)- 

I.e. eigenvalues corresponding to eigenfunctions even or odd with respect to xi 



''A 

h'^h~^ , or as a 2 X 2-system with A„ replaced by I " . 1 and with 2x2 matrix 
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Proof. Proof immediately follows from Bohr-Sommerfeld approximation for A„. □ 

So far propositions 14.51 14.61 were proven for W = 1 but they obviously hold for any 
W > e with the critical values ±1, k* replaced by iVT^/^, k*W^^'^ respectively. Then one can 
easily recover all the results similar to those of subsection 2.3 and thus one can set Ti = e, 
To = Ch~^h\ \ogh\ in inner but not periodic zone Zinn \ ^pcr = {^o < 1^2 — A;*iy2 1 < Cq}. 
These propagation results imply immediately estimate 

|i^t^/^-vxT(^)^(g^n)| < ch"''h~'~'''^\\ogh\ < ch-'^. 

However one can engage the successive approximation method as in subsection 4.1 (with 
p = 1 now); then one needs to consider only the first term of it and following arguments 
of subsection 4.1 one can see easily that \Ft^h^iT-XT{t)'i^{Qipu)\ < Ch^^'j. Thus estimate 
immediately leads to 

Proposition 4.7. For fixed eo contribution of zone Z^^^ \ Z^^r = {eo ^ — k*W^ < Cq} 
to the remainder estimate does not exceed Ch'^'y as fi < eh~^ with small enough constant 
e = e(eo)- 

4.3 Periodic zone. I. Estimates 

Now I need to treat periodic zone Z^^i = {|^2 ~ < cq}. However now I have just one 

equation due to proposition 14.51 and the construction after it. 
We need 

Proposition 4.8. Let ao(^2; K) he a model operator with the potential V = 1 and h < e = 
e(eo). Then 

(4.13) 5|A„,(e2, h) > ei : 16 - ^1 < eo Vn : |A„(6, h)\ < e. 

Proof. Proof follows from Bohr-Sommerfeld approximation and propositions II .41 and II .81 for 

even and odd u respectively. □ 

Now one can apply our standard rescaling technique as in section 3 (but in the simpler 
form) with unperturbed operator A„(/i-D2, ^) + | — ^W{x2) leading to 

Proposition 4.9. Contribution of zone Zp^r = {|6 — k*W^^'^\ < eo} to the remainder esti- 
mate does not exceed Ch~^'y as non-degeneracy condition \{2. 105)~^ is fulfilled and Ch~^^^^ 
otherwise with arbitrarily small exponent S > 0. 
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4.4 Calculations 

Now let us apply our standard successive approximation technique as in section 3 with an 
unperturbed operator SiQ{hD2, h) + ^ — ^W{y2) in the outer and inner non-periodic zones 
and with an unperturbed operator \n{hD2,h) + | — ^W{y2) in the inner periodic zone 
leading exactly to 

Proposition 4.10. Under condition i \4-ll with the remainder estimate equal to Ch~^'y as 
non- degeneracy condition \(2. 105)^ is fulfilled and Ch'~^^^''y otherwise with arbitrarily small 
exponent 6' > asymptotics 

(4.14) / e{x,x,0)ip{x2) dx = {27ihy^ / no{x2,^2,0)ip{x2) dx2d^2- 



Since S^"^ = S^"^ mod 0(1) under condition ()4.1|) this proposition finalize the proof 
of Theorem 10.11 under condition ()4.1|) . 

4.5 Ultra-strong magnetic field 

Now we are left with the case 

(4.15) eh-" <fi< Coh"" 

when h is not a small parameter anymore and when all the eigenvalues of ao(^2, K) which 
are less than C are separated by a small constant. 

Then one can reduce the analysis to those of 1-dimensional scalar operators A„(x2, hD2, h) 
and the problem is that we don't know how A„(^2; ^) depends on ^2- 

So far I can prove unconditionally only 

Proposition 4.11. Under condition 

(4.16) \Xn{^2,h) + ^~^W\ + \d,,W\>e, Vn,e2,X2 



asymptotics [j.l^] holds modulo 0(1). 



Note that for even v perturbations Bnn{x2i^2) are 0{h'^) since terms of the order h 
could appear only from multiplication by Xi in the first power but due to the parity of 
eigenfunctions these terms are 0. For odd v we can conclude that i3ri„(x2,^2) = Oiji^) as 
^2 = but I want to make conjecture 

(4.17) S„n(a^2,6) = 0{h^) as = <9,.,A„ = 
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since I don't know if there are other critical points but 0. 



Then according to |Ivrlj . section 4.4 the remainder estimate should be 0(1) under 
conditions 



(4.18) 



E 

l<k<l 



with / = 2 and nondegeneracy condition (j2.1(J3|) . Further, note that condition |(4.18);| is 
fulfilled with some / due to the analyticity of A„(^2, ^) with respect to ^2 and the fact that 
An(^2,^) — ^ 00 as 1^2! 00. Then according to jlvrlj . section 4.4 remainder estimate 
0{h~^ ) in the general case is guaranteed but to avoid contribution of junior symbols into 
the principal part one should make additional assumptions. 



Proposition 4.12. (i) Under conditions d^.i^ , ( |//.i7D , (4.18)2 and <\2. 1 U'^ asymptotics 
( 4-14 ) holds with the remainder estimate 0(1). 

(a) Under conditions { 4-151 , ( ^. j7 ) and (4.18)2 asymptotics { 4-14 ) holds with the remainder 
estimate 0{h~^). 



However analysis of to jivrlj . section 4.4 shows that the the remainder estimate must 
be 0{h~^) under [(4. 18);| with any / and non- degeneracy condition [(2.105)"^ with any m but 
there would be indefinite correction of magnitude 0(/i^/^"'"^/™'~^). Further, under condition 
(2.105)2 remainder estimate is 0(1). Therefore we arrive to 



Proposition 4.13. (i) Under conditions Ul^ , |(4.18)f| with I > 3 and (2.105)2 

asymptotics { 4-14 ) holds modulo 0(1); 

(a) Under conditions i \4-15^ , i \4-ni , ^4:.18)i\ with I > 3 and \{2. 105) m\ with m > 3 asymptotics 
dCTl ) holds modulo ojW+^/^^ . 



5 Generalizations 
5.1 Vanishing V 

We can also generalize all our results to the case of V vanishing: 
Theorem 5.1. Theorem W. 1\ remains true if condition (1(7. gjl is replaced by 

(5.1) \d.,V\>eo. 
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Proof. Consider ^-admissible partition with £ = e\V\. 

(i) First, consider elements on which £ < ei|xi|. Then rescaling x ^—>- x/i and diving by £ 
we arrive to the situation of non- vanishing magnetic field fics = A^Ti"^-^^''^ with Fcs = F'y~'^ 
and Ves = V£~^ with a parameter 7 x |xi|. 

Further, after rescaling Ves/Fes = (V/F) ■ (Y/^) ^ 1 and also £| Wcff/i^cffl x 1- 
To apply "non-vanishing" theory with hes = h£~^^'^ and (Xes one must assume that 
^cff < 1, fJ'cS > 1 or equivalently 

(5.2) £ > Comax(/ii, {fiY~^y^) = ^■ 

Then the contribution of such element to remainder estimate does not exceed Cfi~gh~g = 

After summation over partition with fixed 7 one gets due to condition (|5.ip the same 
expression but with £ = 7 i.e. Cfi~^h~^'y^~'^. One must assume that condition (|5.2|) is 
fulfilled for this top £, which is equivalent to 

(5.3) 7 > Comax(/ii,/i-2/(2-i)) = L. 

Then summation over a;i-partition satisfying this condition results in fi'^h'^^L'^'" which as 
one can check easily does not exceed Cfi~^^'^h~^. 

(ii) One needs also to consider elements satisfying ()5.3j) but not ()5.2j) . Let us unify such 
elements, so in fact I redefine £ setting it £ = e|l^| -|- £. Then on elements with £ x £ either 
/ieff X 1 or /igfT ^ 1- In the former case the condition of a potential being disjoint from 
is not needed while in the latter case element in question is forbidden as fies > C*i and 
the condition of a potential being disjoint from is not needed again. The contribution to 
the remainder estimate of of each of these elements is either Ch^^^ (as yUefr ^ 1) or C/i^ 
(as /icflf X 1) and in any case it does not exceed Cn'^hc^ again. Repeating summation 
procedure of (a) one can see that the total contribution of such elements does not exceed 

So one needs to consider two other types of elements. 

(iii) Consider first elements on which 

(5.4) \xi\<C£, £ = e\W\. 

Then rescaling x ^ x/ £ and dividing by £ brings us to the frames of our main theory with 
/ieff = jJi£''~^^'^, = M"'^/"^ and one must assume that 

(5.5) Comax(/ii,/i"2/(2--i)). 
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Then contribution of each element to the remainder estimate does not exceed Cfi'^ h~Q = 
C ix~^^'^h~^i^'^^^^^'^ and the sum over partition results in Cfj.~^^'^h~^. The exceptional ele- 
ments will be covered in the next paragraph. 

(iv) So, we are left with the analysis of the zone {|xi| + \W\ < L} with L defined by 
fl5.3p . Due to condition ()5.1|) with no loss of the generality this zone could be replaced 
by < L}. Rescaling x ^ x/L and dividing by L we arrive to the same situation as 
we started but with = hL~^^'^, fi^s = iiL^~^ instead of h^jj, i.e. with either /lofr < 1, 

/ieflF X 1 or with /icff X 1, /Ucff > 1- 

In the former case I just refer to the classical theory: contribution of this zone to the 
remainder estimate does not exceed 0{h~^ = Ch~^L^/'^ < Ch~^ jjL~^/'^ . In the latter case 
contribution of this zone to the whole asymptotics does not exceed C^~^ as this zone is 
"forbidden"; one can prove it easily going to the axillary space //^(M^-J. □ 

5.2 Remark about higher dimensions 

3-dimensional case is rather boring because the best remainder estimate if magnetic field 
does not degenerate is 0{h~'^) and applying the standard rescaling arguments one will 
arrive to 0{h~^ j dx) and in the generic case 7 is the distance to the finite set, thus 
resulting in 0{h~'^) again (as /i < h~^)). As fi > and \F\ x 7 the principal part will 
be of magnitude h~^'yf and the remainder estimate of magnitude h~'^'jf with 71 = {fih)~^ 
as ^ fJ^ < Ch'"^. One can expect the similar results in higher dimensions too. 

Even-dimensional case (I will treat it in one of forthcoming papers) is much more in- 
teresting but still less than the special 2-dimensional case. First of all as it follows from 
Martinet |Maj in the generic case as d > 2 magnetic matrix F never vanishes and thus 
Tr^ g~^F > (where Tr"*" g~^F is the sum of positive eigenvalues of ig~^F and is a met- 
rics matrix); therefore unless we replace V hy V + Tr"*" F the principal part will be just 
and the remainder estimate 0(/i~*) as yU > Ch~^ thus cutting us from the case in which 
one expects the periodic zone to be important. 

Canonical form of F in the generic case is known due to |Maj and Roussarie; first of 
all if /i is the smallest positive eigenvalue then /i x xi in the appropriate coordinate 
system. Further, in dimension 4 F has no more than 2 degenerating eigenvalues while 2 
others are disjoint from and thus 4-dimensional case seems to be nothing but the lift-off 
of 2-dimensional (as /i < Ch~^). It is not that simple in particular because there appears 
a very special submanifold of codimension 3, but likely the above assertion is correct. 

In dimension d>Q more than 2 eigenvalues can degenerate but the degeneration hap- 
pens on the manifold of codimension 6 and it seems to be too thin to have any effect. It 
does not mean that the proofs are going to be easy but that one should not expect any 
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surprises. 
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